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A minimal model for magnetic reconnection and, generally, low-frequency dynamics in low-beta 
plasmas is proposed. The model combines analytical and computational simplicity with physical 
realizability: it is a rigorous limit of gyrokinetics for plasma beta of order the electron-ion mass ratio. 
The model contains collisions and can be used both in the collisional and collisionless reconnection 
regimes. It includes gyrokinetic ions (not assumed cold) and allows for the topological rearrangement 
of the magnetic field lines by either resistivity or electron inertia, whichever predominates. The 
two-fluid dynamics are coupled to electron kinetics — electrons are not assumed isothermal and 
are described by a reduced drift-kinetic equation. The model therefore allows for irreversibility 
and conversion of magnetic energy into electron heat via parallel phase mixing in velocity space. 
An analysis of the exchanges between various forms of free energy and its conversion into electron 
heat is provided. It is shown how all relevant linear waves and regimes of the tearing instability 
(collisionless, semicollisional and fully resistive) are recovered in various limits of our model. An 
efficient way to simulate our equations numerically is proposed, via the Hermite representation of 
the velocity space. It is shown that small scales in velocity space will form, giving rise to a shallow 
Hermite-space spectrum, whence it is inferred that, for steady-state or sufficiently slow dynamics, 
the electron heating rate will remain finite in the limit of vanishing collisionality. 



I. INTRODUCTION 

Multiscale nonlinear plasma phenomena pose 
some of the most tantalizing and intellectually chal- 
lenging theoretical problems in the field. Roughly 
speaking, this is because they tend to involve cou- 
pling between fluid (large) and kinetic (small) scales 
and so strain both our physical intuition (which in 
most cases is anchored in the fluid description) and 
analytical stamina (which tends to break down when 
full Vlasov-Maxwell kinetic theory appears to be the 
only rigorous recourse). A particularly clear exam- 
ple of this situation is the development of the (non- 
linear) kinetic theory of magnetic reconnection. 

Magnetic reconnection, the unfreezing of magnetic 
flux in nonideal conducting plasmas is one of 
those fundamental plasma physical processes under- 
standing which turns out to be an essential theoreti- 
cal building block in a vast range of laboratory, space 
and astrophysical applications: the sawtooth crash 
in tokamaks^ solar flares^ 5 - the magnetosheath 
and magnetotail of the Earth^r— various types of 
magnetohydrodynamic and plasma turbulence^^ 
particle acceleration mechanisms 1 1 1 1 1 2 as well as a 
plethora of more exotic reconnection-related phe- 
nomena believed to occur in extreme astrophysical 
environments^ — solutions to these problems turn 
out to depend crucially on whether one thinks mag- 
netic reconnection is fast or slow, steady or bursty 
and impulsive^ collisional or collisionless ) 15 ! 16 es- 
sentially two- or three-dimensional, mediated by 
resistivity ; 17 ' 18 dispersive waves^ instabilities of the 



current sheets^ — turbulence^ 6 - — or various com- 
binations of these. 

There appears to be a broad theoretical consen- 
sus that in many natural plasmas, reconnection can- 
not be fully understood in a purely magnetohydro- 
dynamic (MHD) setting (the same is true for most 
other nonlinear multiscale plasma phenomena). In 
recent years, a dramatic increase in the raw comput- 
ing power that could be brought to bear on these 
problems has fuelled an intense effort — and some 
spectacular empirical progress — in understand- 
ing kinetic reconnection using PIC simulations ^2r— 
both in two and, more recently, three dimensions. 
One limitation of this approach is the lack of a sim- 
ple analytical framework that can help in the inter- 
pretation of numerical results. Ever since the pure 
MHD approach was recognized as insufficient, and 
even before that, attempts have continued to de- 
velop various minimal models, usually of the two- 
fluid kindr 6 - - — with the aim of capturing the "es- 
sential physical ingredients," reflecting at any given 
time the evolving understanding of what those were. 
Finding such minimal models has been and is likely 
to remain necessary also because, despite the rapid 
expansion of brute-force simulations, the resolu- 
tion available, however impressive, is, in fact, never 
enough: the problem has at least three well sep- 
arated spatial scales (global fluid, ion kinetic and 
electron kinetic), as well as possibly a similar level 
of complexity in the kinetic phase (velocity) space, 
and is likely to require three dimensions in an essen- 
tial way. Furthermore, the current fully kinetic PIC 
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simulations 3 -^ have as yet to produce a clear picture 
of reconnection in low-beta plasmas embedded in a 
strong guide field, while gyrokinetic simulations of 
magnetic reconnection 4 ^— are in their infancy and 
are afflicted both by high numerical resolution re- 
quirements and lack of clarity on the theoretically 
expected outcomes (the gyrokinetic theory 4 ^ in its 
general form is not much less analytically difficult 
that the full Vlasov-Maxwell kinetics, although nu- 
merically it is, of course, a radical simplification). 

There is a good reason to believe that a physically 
realizable model set of equations for weakly colli- 
sional reconnection, or nonlinear plasma dynamics 
generally, however simplified, cannot in general be 
fluid. Like most nonlinear phenomena, magnetic re- 
connection can and, indeed, should be thought about 
in terms of energy conversion, namely, conversion 
of magnetic energy associated with the reconnecting 
large-scale configuration into other forms of plasma 
energy and, ultimately, particle heat. It is when the 
energy is dissipated into heat that this conversion 
becomes irreversible. However, when collisions are 
weak, Ohmic (resistive) heating is not an important 
process, so most purely fluid models of collisionless 
reconnection turn out to be Hamiltonian; 36 ' 37 i 40 ' 41 
reconnection converts magnetic energy into other 
forms of "fluid" energy (e.g., kinetic energy of the 
mean electron flows), the entropy of the system does 
not increase, everything is in principle reversible, 
and there is no heating. This is somewhat similar 
to inviscid dynamics in a neutral fluid. In fact, just 
like in neutral fluids, nonlinear dynamics in plasmas 
generically trigger formation of small scales, so even 
small dissipation coefficients become non-negligible 
because they multiply large gradients (which is why 
dynamics in fluids with small viscosity are not ev- 
erywhere inviscid, boundary layers form). In weakly 
collisional plasmas, the small-scale structure arises 
in phase space (i.e., both in position and particle ve- 
locities) via linear and nonlinear phase mixing pro- 
cesses (see Sec. IV Dl) . Once large velocity-space gra- 
dients are present, even weak collisions are sufficient 
to dissipate energy and cause heating. Formally, 
this can be understood by tracing the evolution and 
flows of free energy — the quadratic invariant whose 
transfer to small scales and eventual thermalization 
is the central process in the non-equilibrium ther- 
modynamics of kinetic plasmas (in the context of 
kinetic an d gy rokinetic turbulence, this is explained 
in Refs . [4811491 where further references are provided; 
in application to magnetic reconnection, we expand 
on this topic in some detail in Sec. lIV[) . Thus, we be- 
lieve there is a need for a minimal model of weakly 
collisional reconnection that is not fully conserva- 
tive, because exact conservation properties impose 
constraints on the phase space of the syste m 50 ' 51 



the breaking of which is not just non-negligible but 
is in fact likely to be physically essential in real 
plasmasi 1 ^ While numerical evidence is perhaps not 
fully conclusive on this subject, a few recent studies 
have, explicitly or implicitly, stressed the importance 
of electron heating in the nonlinear regime of kinetic 
reconnection i 35 i 45 ' 52 

In this paper, we take the view that a minimal 
model that is as fluid-like as possible is clearly desir- 
able, but ad hoc closure approximations are danger- 
ous even if they appear to be physically motivated. 
We would like therefore to have a model that com- 
bines analytical and computational simplicity with 
physical realizability, i.e., constitutes a rigorously 
correct approximation of the kinetic system in some 
well defined physical limit. It turns out such a model 
can indeed be constructed and is a simple general- 
ization of an existing two-fluid modelj^l although 
it is not a fluid model in that it does retain kinetic 
electrons: the equations are the continuity equation, 
the "gyrokinetic Poisson equation" for the ions, the 
generalized Ohm's law and a version of the drift- 
kinetic equation for electrons; the latter is coupled 
to the fluid variables via the gradient of parallel elec- 
tron temperature (energetically, via work done by 
the parallel electron pressure gradient). The formal 
limit in which this model is derived is a combina- 
tion of the gyrokinetic regime (strong guide field, 
frequencies below ion cyclotron, strong anisotropy 
fcii <C kj_, small-amplitude fluctuations) with a low- 
beta expansion (plasma beta of order electron-ion 
mass ratio) — this is explained in detail in Sec. Ill Bl 
A very close precursor of our approach in the existing 
literature is the fluid-kinetic model proposed by de 
Blank ; 53 ' 54 which couples a two-fluid model to a sim- 
plified electron kinetic equation, although he does 
not give a rigorous asymptotic ordering under which 
his model holds and also restricts his attention to the 
exactly 2D, exactly collisionless case, which means 
that his model is Hamiltonian and has an infinite 
number of Lagrangian conserved quantitie s 55 ' 56 (see 
Appendix [A} . 

Let us discuss the basic physical ingredients that 
are retained in our approach. Despite the simplicity 
of our equations (summarized in Sec. IIII E[) . we be- 
lieve that very little is lost of what we consider the 
sine qua non of kinetic reconnection. 

Three scales. The problem is fundamentally 
three-scale: the model has to allow for a reconnect- 
ing configuration on fluid scales, dispersive effects 
at the ion scales and the flux unfreezing by a colli- 
sionless mechanism associated with electron micro- 
physics and so residing at electron scales. In our 
case, the ion scales are the ion sound and Larmor 
radii (the ions are gyrokinetic and not assumed cold; 
see Sec. IIIIC|) . The flux unfreezing is effected by 
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electron inertia — see Sec. IIII Bt we will sacrifice 
the electron-Larmor-radius effects. 

Collisions. It is desirable for a good model to 
contain a smooth transition from collisionless to 
semicollisional to fully collisional (resistive MHD) 
regime — both because this provides a way to 
benchmark against situations that are better under- 
stood and because theoretically it is possible that 
many natural systems teeter at the boundary be- 
tween the collisionless and collisional regimes . 15 ' 16 
The transition between the two has been the focus 
of several recent studies, both experimental 5 -^ and 
theoretical) 34 ' 35 ' 46 ' 58 but remains poorly understood. 
In our model, electron-ion collisions are retained (see 
Sec. IIII D]) and so both Ohmic resistivity and (paral- 
lel) electron heat conduction (in the semicollisional 
limit; see Sec. IV C[) are recoverable (but not ion or 
electron viscosity, a limitation that will be discussed 
further in Sec. IVII). 



Free- energy flows, Landau damping and electron 
heating. Finally, as we explained above, the model 
provides an electron heating channel, operative 
even with very weak collisions, in the form of 
the coupling of the generalized Ohm's law to the 
electron kinetic equation with a collision operator 
retained. This means that the electron Landau 
damping (whose essential byproduct is parallel 
phase mixing) is included and free energy can be 
converted from various fluid forms into the electron 
entropy and thence to heat (see Sec. IIV Bp . Note 
that neither ion Landau damping nor the nonlinear 
perpendicular phase mixing of either specie s 48 ' 49 ' 59 
survive in our model. 

The rest of the paper is organized as follows. In 
Sec. UH we review the gyrokinetic system, which 
is our starting point (Sec. [TT~^|) . and introduce the 
set of ordering assumptions that encode the phys- 
ical limit in which our equations hold rigorously 
(Sec. IIIB|) . In Sec. IIIII we derive the equations 
themselves (they are summarized in Sec. IIII Ep . 
In Sec. IIVI we work out the energetics of these 
equations: the various forms the free energy takes 
(Sec. IIV A[) . how it is exchanged between fields 
(Sec. IIV Bj) and what that implies about irreversibil- 
ity, thermalization (dissipation), and electron heat- 
ing (Sec. IIV Cp . In Sec. [Vj we introduce a spec- 
tral representation of the electron kinetics in terms 
of Hermite polynomials (using a modified Lenard- 
Bernstein operator for electron collisions, introduced 
in Sec. IIII D[) . This provides what appears to be 
both a remarkably simple computational approach 
and an intuitively appealing physical interpreta- 
tion of velocity- space dynamics and electron heat- 
ing as a cascade in Hermite space (Sec. IV Dp . We 
are led to the conclusion that the electron heat- 



ing rate should remain finite in the limit of pos- 
itive but vanishing electron collision frequency - 
except for fast-growing solutions like the tearing 
mode (Sec. IV Ep . We also derive the semicolli- 
sional limit of our equations 6 ^ in Sec. IV CI Sec- 
tion IVII contains the concluding discussion. The 
paper is supplemented with two technical appen- 
dices: on two-dimensional invariants of our sys- 
tem (Sec. and on the linear theory (gyroki- 
netic electron-Landau-damped Alfven waves^i colli- 
sional (resistive MHD)(£2r— semicollisional 6 ^— and 
collisionless 6 ^ 67 ' 71 " — tearing modes) . 

We have adopted a rather gradual, step-by-step 
approach to the derivation of all results. An impa- 
tient reader, or one already familiar with most of 
the (standard) analytical machinery deployed here, 
can gain a basic idea of the main results and con- 
clusions by looking at Sec. BlTEl Sec. HVCl Sec. EH 
and Sec. EU 



II. PRELIMINARIES 

A. Gyrokinetics in a slab 

The starting point for our derivation is the gyroki- 
netic description of magnetized plasma*^ which is 
appropriate for low- frequency (tu -C Q s ) anisotropic 
(fc|| <C fe_i_) fluctuations in the presence of a mean 
magnetic field. The simplest case, which is all we 
will require here, is that of a static uniform equilib- 
rium with zero electric field and a constant straight 
magnetic field Bo = BqZ, whose direction defines the 
z axis. The subscripts _L and || will always refer to 
directions with respect to this equilibrium ("guide") 
magnetic field. The gyrokinetic equations are de- 
rived by performing an expansion of the Vlasov- 
Landau equation in the small parameter e = 
and averaging out the particle Larmor motion, i.e., 
all frequencies greater than the cyclotron frequency 
Sl s = q s Bf)/m s c, where s = i, e is the species index, 
q s is particle charge, m s particle mass, and c the 
speed of light. For a detailed derivation of homoge- 
neous gyrokinetics in a slab, the reader may consult 
Ref. HI 

Let us summarize the resulting equations. The 
distribution function up to first order in e is 

f a (r, v, t) = F 0s - ^P^-Fo a + h s (R Sl v ± ,vu , t), 



Tos 



where the zeroth-order distribution 



Fos(v) = 



n 0s 



3/2 



exp 



"ths 



(1) 



(2) 
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is a Maxwellian with uniform density no s and tem- 
perature Tq s , w t hs = (ZTos/ms) 1 / 2 is the thermal 
speed, q s (p/To s = 0(e) is the Boltzmann response 
containing the electrostatic potential ip, and h s is 
the gyrocenter distribution function, also 0(e). Spa- 
tially, h s is defined as a function of the the gyrocen- 
ter, or guiding-center, coordinate 



R s = r 



a, 



(3) 



and satisfies the gyrokinetic equation 

dh s dh s c , , 



qsF 0s d(x) Rs fdhs 

t 0s at \ at 



(4) 



B. Low-beta ordering 

We would like to derive a minimal model suitable 
for an analytical description of magnetic reconnec- 
tion in the presence of a mean field, i.e., reconnec- 
tion of antiparallel perturbations SB± = — z x V>1||. 
Since we wish to have a model that describes a real 
physical situation, we cannot resort to writing ad 
hoc fluid equations. Instead, our model will take the 
form of an asymptotic expansion of the gyrokinetic 
equations under an appropriate physically motivated 
ordering of all spatial and time scales and of the per- 
turbation amplitudes. In devising our ordering, we 
are guided by what is known or expected about the 
physical effects that are essential in any description 
of a kinetic reconnection process. 



where x(r, v, t) = ip — v ■ A/c is the gyrokinetic po- 
tential (containing the information about the elec- 
tromagnetic field in the form of the scalar potential 
ip and vector potential A of the perturbed magnetic 
field, SB = V x A), 

is the Poisson brackets, (dh s /dt) c is the (gyroaver- 
aged) collision operator, and 



(X(r,v,t)) 



1 

2^ 



Vj_ x z 



-, V 



(6) 

is the average of \ at constant R s over the gyroan- 
gles "d. In Fourier space, this gyroaveraging opera- 
tion takes a simple mathematical form in terms of 
Bessel functions Jo an d Ji , so the Fourier transform 
of (x(r, v, t)) R with respect to R,, can be written 
as follows 



(x) Rs ,u = J o(a s ) ipk 



vuAu 



T 0s 2v 2 ± Ji(o.) (5B||k 



q s v. 



ths 



(7) 



where a s = k±vj_/tt s , and tp^, An^ and 6B^ are 
Fourier transforms (with respect to r) of the scalar 
potential, parallel component of the vector poten- 
tial and parallel component of the perturbed mag- 
netic field, respectively. These fields are determined 
via Maxwell's equations, namely, the quasineutrality 
and Ampere's law, where particle densities and cur- 
rents are calculated from the gyrocenter distribution 
h s . These equations will be introduced in Sections 
MI Al and lHI CI where we will need them to compute 
electron flow velocity and density perturbation. 



1. Spatial scales 

Firstly, in the presence of a strong guide field, elec- 
tron inertia is expected to be a key mechanism for 
breaking the magnetic field lines (flux unfreezing) 
in collisionless or weakly collisional plasmai 62 ' 77 " — 
Thus, we order 



k±d e ~ 1, 



(8) 



where d e 



c/ujp e 



P&l\fWe is the electron in- 



Pe = Vthe/tte the 

8Trn 0e T e/B 2 the 



ertial scale and, to fix standard notation, 
(47rno e e 2 /m e ) 1 / 2 is the electron plasma frequency, 
e = \q e \ the elementary charge 
electron Larmor radius, and /3 e 
electron beta. 

Secondly, a key feature of kinetic reconnection is 
a double layer resulting from the decoupling of the 
electrons from the ions at the ion sound scaled To 
retain this effect, we order 



1, 



(9) 



where p s = pi ' Z/2t is the ion sound radius (sub- 
script s for "sound" not to be confused with the 
species index!), pi = vthi/^i is the ion Larmor ra- 
dius, Z = qi/e and r = T i/Toe- We will consider 
the temperature ratio to be order unity, so Eq. (|9|) 
immediately implies 



1. 



1/2 



1. 



(10) 



Thus, we retain the ion FLR along with the ion 
sound scale. However, a further simplifying option 
remains open: assuming cold ions, t « 1, we can 
eliminate the ion FLR effects. It is with a view to 
this possibility that we will carry the r dependence 
in all our orderings discussed below. 
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In order for Eqs. © and © to be consistent, we 
must have, within our ordering, d e ~ p s (this does 
not mean that these scales must be similar, just that 
we are not hard-wiring into our model a disparity 
between them). This implies 



de 
Ps 



<2Z 



1/2 



To achieve this, we order 



A 



Znip 



m. 



< 1. 



(11) 



(12) 



Thus, we are restricting our consideration to low- 
beta plasmas and allowing both species to have fi- 
nite temperature. The ordering (|12|) is appropri- 
ate, for example, for the solar coronal and low-beta 
laboratory experiments such as the LArge Plasma 
Device at UCLA^i In modern tokamaks, values of 
Pe ~ 10 -3 can occur in the edge pedestal region in 
the H-mode regime^ 

Note that in the gyrokinetic approximation, /3 e is 
considered order unity with respect to the ordering 
of all quantities in powers of e = k\\/k±; we will 
treat our low-beta expansion as subsidiary to the 
gyrokinetic e expansion. 

Now, using Eq. ijHJ), we conclude 



k±p e 



< 1, 



(13) 



which will allow us to neglect the electron FLR ef- 
fects and derive an "almost fluid" set of equations for 
the electrons. 



2. Time scales and perturbation amplitudes 

Far from the reconnecting region, the plasma mass 
flow is ordered with the E x B drift velocity. Thus, 
we order the fundamental time scale on which we 
allow our fields to vary in such a way that the char- 
acteristic frequency is that associated with the E x B 
velocity uj_: 



kl 



B~ 



(14) 



Note that the gyrokinetic approximation requires 

Since it is an essential feature of our model to 
take into account electron kinetics, we must allow 
the parallel streaming frequency of the electrons to 
be the same order as the rate at which our fields 
vary, namely 



The requirement that Eqs. (fT4"|) and (fT5j) should 
be consistent with each other imposes an ordering 
on the size of the scalar potential: 



'It 



1 



eip 

T e k± kj_p e 



(16) 



where we used Eq. (fTU|) . Note that the appearance of 
the gyrokinetic expansion parameter e in the order- 
ing of the perturbation amplitudes confirms that the 
use of the gyrokinetic approximation is appropriate 
in the physical circumstances we are considering. 

We further require that the density perturbations 
are of the same order as the electrostatic perturba- 
tions given by Eq. (IK 



Sn e 
noe 



Z eip 
t T 0e 



Z 



(17) 



(the factor of Z/t will emerge in Sec. MI CI and is 
kept for book-keeping purposes) . Physically this fol- 
lows from the requirement that the physics associ- 
ated with the ion sound scale is retained [Eq. ([9])]. 



3. Alfvenic perturbations 

We will now order the magnetic perturbations. 
Let us observe that the ordering (fl"2l) implies that the 
electron thermal speed is comparable to the Alfven 
S /\/47rm l n 0l : 



speed va 
(quasineutrality) , we have 



«the 

va 



/3 e rrij 

Z TO, 



since noi 



1/2 



Therefore. 



■\\VA, 



n 0e /Z 



(18) 



(19) 



i.e., Alfven waves can propagate along the guide field 
with the same characteristic frequency as the elec- 
trons stream and plasma flows. Note that in view of 
Eqs. <j9j) and (fT0|) . this ordering holds both for the 
magnetohydrodynamic Alfven waves (oj = k\\VA), 
and for the kinetic Alfven waves (lu ~ k^VAk±p s ). 

Stipulating that Alfvenic perturbations should be 
accommodated by our ordering (which is equivalent 
to demanding that the Lorentz force is nonnegligible 
— an essential ingredient in a reconnecting system) , 
we deduce the ordering for the perpendicular per- 
turbed magnetic field <5Bj_ = — z x V±A»: 



5B ± u ± (Z 
e k xPi _ 

Bq v a \t 



1/2 



(20) 



k\\V t he 



(15) where we have used Eqs. (|HT|) and (JT 



G 



Note that, since e ~ k^/k±, this implies k±SB± ~ 
fc||i?o, i-e., the spatial variation of all quantities along 
the exact magnetic field contains comparable con- 
tributions from their variation along both the mean 
and perturbed fields — a general property of gy- 
rokinetic perturbations. This is just what is needed 
for reconnection problems with a guide field because 
Alfvenic dynamics with respect to the perturbed 
field must be allowed. 

Finally, from the perpendicular component of Am- 
pere's law (see, e.g., Eq. 120 of Ref . 1491) . 



6B\\ 



eip 



(21) 



where Eq. (|T6"| has been used. This will cause the 
parallel perturbations of the magnetic field (i.e., per- 
turbations of the field strength) to fall out of our 
final set of equations. Indeed, it is a well known 
fact that in the low-beta ordering such perturbations 
tend to be negligible. 

This completes the ordering of the perturbation 
amplitudes. 



where we have used Eqs. (fT2"j) and (f2"2"T) . Thus, we 
are effectively assuming ions to be collisionless (no 
ion viscosity). 

Finally, it is perhaps useful to note the ordering 
of the particle mean free path: using Eqs. (fT5)) and 
(J22]), we have 



&|| A m fpe 



mfpi 



k\\vth, 



(26) 



III. DERIVATION OF THE EQUATIONS 

We are now ready to apply our ordering and de- 
rive from the gyrokinetic system (Sec. Ill A[) a re- 
duced system of equations describing gyrokinetic 
ions and drift-kinetic electrons in a low-beta weakly 
collisional plasma. An impatient reader can skip to 

sec. urn 



Electrons 



4- Resistivity and collisions 

While our main focus is on collisionless reconnec- 
tion, we would like to make contact with the col- 
lisional limit, in which the magnetic flux unfreez- 
ing and magnetic energy release is accomplished by 
Ohmic resistivity^ 3 - We can retain resistivity by or- 
dering the electron-ion (and consequently electron- 
electron) collision frequency as comparable to the 
characteristic frequency of all the other processes 
that we are taking into account: 



Zv e 



(22) 



Since the Ohmic magnetic diffusivity (often colloqui- 
ally called resistivity) is r\ ~ Veid 2 , the above order- 
ing means that the diffusive effects are retained in 
our ordering: indeed, using Eqs. (jHJ and (|22jl. 



rjk 



■k 2 d 2 



(23) 



The resistivity can later be neglected in a sub- 
sidiary collisionless expansion, but we consider it 
useful to have a model in which a smooth transi- 
tion from collisional to collisionless regime is possi- 
ble (cf. Refs. [34.35 46 57] and ideas on the marginal 
collisional-collisionless reconnection regim o ' 1 ) . 

Our ordering of the electron collision immediately 
implies an ordering of the ion collisions <£i 



mi 



1/2 



3/2 



(24) 
(25) 



In Eq. ([5]) for electrons (s = e), we retain only 
the lowest order in the expansion with respect to 
f3 e . In order to do this, we note first that the Bessel 
functions in the expression for (x) R [Eq. ([7])] can be 
expanded in small argument because a e ~ k±p e ~ 
VAT< 1 [Eq. (USD]: 

j o{ae)c ,^lM = 1 + (a 2 e ). (27) 
a e 

Using Eqs. (jTr?)) and (|2"0|) , we find that the ip and An 
terms in Eq. ((7|) are the same order, 



vii A11 v the B 5B A 



T 0e 



ck 1 



Bn 



e k±p e 



(28) 



while, according to Eq. (f2"Tj) . the 6B\\ term is one 
order of /3 e smaller. Thus, 



(X) 



R, 



1 A, 



[1 + 0(/3 e )] . (29) 



This allows us to write Eq. Q as follows, up to 
corrections of order 0((3 e ), 



dh e 
~~dt 



b-V/i P = 



eFpe d 
T 0e dt 



v\\A\ 



dhe 
dt 



(30) 

where we have introduced a "convective" time deriva- 
tive incorporating the E x B motion and the parallel 
spatial derivative along the perturbed field line: 

dh e dh e c , , 

-# = -m + B- ^ h ^ (31) 

dh 1 

b-V/ le = -Qf-B- { A \\M- (32) 
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Note that all terms in Eq. ([50)1 are comparable under 
ordering because of the assumptions we made about 
the time scales for the E x B flows [Eq. (fT4[) [. elec- 
tron streaming [Eq. (fT5"[) [. collisions [Eq. (J22)] and 
perturbation amplitudes [Eqs. ([To]) and ([2"0"[) [. 

It is now convenient formally to split the electron 
distribution function in such a way as to separate the 
inhomogeneous solution arising from the first term 
on the right-hand side of Eq. ([3"0"[) as well as the 
density and parallel velocity moments: 



^ + — + -72 }F0e + 9e, (33) 



J the 



where by definition 

5n e 1 
n-oe n 0c 
1 



d 3 v6f e 



= / d VVnSfe, 



(34) 
(35) 



and 6f e = h e +etpFo e /To e is the total perturbed elec- 
tron distribution function [see Eq. (UJ]. Since the 
first term in Eq. (j33[) cancels the Boltzmann part in 
(5/ e , the above definitions of the density and paral- 
lel electron flow velocity are consistent provided we 
demand that 



' v\\ ' 



(36) 



i.e., the "reduced" electron distribution function g e 
contains all higher moments of the electron distribu- 
tion function, but not density or parallel flow. 

The perturbed density will be calculated in 
Sec. MI CI from the quasineutrality condition. The 
parallel electron velocity ui\ e is related to An and to 
the parallel ion velocity u\u via the parallel compo- 
nent of Ampere's law: 

A<rr 

^•(ViXiBi) = -vi^n = — j|| 

c 

4iren 0e , . 

= \U\\i - 1i|| e J . (37) 

This can be easily manipulated into 

«|| e = — d 2 e W\A\\ +u h . (38) 
cm e 



We now substitute Eq. ([33]) into Eq. ([30]) and take 
moments of the resulting equation. The zeroth mo- 
ment is, after using Eq. 



d Sn e 
dt n 0e 



b-Vu 



II' 



b-V— ^Vi^|,. (39) 



The first v\\ moment of Eq. ([3"U[) is, again using 
Eq. (|38p and dividing through by eno e /cm e , 



d 1 A 
dt^ 

d<p 



dlVUw) 

cT 0e 



dz 
cm e 
enoe 



ST 



r/'S 



dh e 
~dt 



To, 
crrir du 



II' 



dt 



(40) 



where the parallel electron temperature perturba- 
tion has been 
moment of g e 



tion has been introduced as a shorthand for the ujj 



ST 



Toe 



1 

riQe 



2v\ 



(41) 



J the 



Finally, an equation for g e is obtained from 
Eq. ([301) by subtracting from it Eq. J3H]) multiplied 
by Foe and substituting dA\\/dt calculated using 
Eq. (|4"0"[) . After a few lines of straightforward al- 
gebra, this gives 




C[9e 



dy±A {] 



(42) 



C[g e 



where the collisional terms have been assembled to- 
gether: by definition, 

'dh e \ 2v\\F 0e f, 3ll (dh e 

(43) 

It is easy to verify that Eq. (|4"!2)) respects the con- 
straint on g e given by Eq. ([3T)[) . We stress that 
Eq. (|42p is not homogeneous ing e , so g e cannot be 
consistently neglected (cf . Ref . [371 ) . 

B. Flux conservation 

Equation (|4T))) can easily be rearranged into the 
following form 



dt 
where 



d A d 

--cb.V^ + r ? ViA||+-dM^lh ( 44 ) 



ip = if ■ 



T 0e ( 5n e 8T\\ e 



n 0e Ti 



(45) 



Recalling that <5B^ = — z x Vj_Aii and B = Bqz + 
5Bj_, we immediately infer from Eq. (|4"4"[) that 

<9B d 

— = V x (u cff x B) + rjVlB -zxVi- d 2 e V\A\\ , 

(46) 
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where u e ff = z x V±clp/Bq is the effective veloc- 
ity, into which, as the above equation demonstrates, 
the field lines are frozen except for the resistive and 
electron-inertia effects. This is an application to our 
equations of the more general flux-conservation ar- 
gument due to Cowley^ 

An important conclusion is that the kinetic elec- 
tron effects, which enter via 5Tji e determined from 
Eqs. (|41]) and (|32]). do not unfreeze flux. 



Ions 



Equations (f3T)H4"2")) are four equations that involve 



ST\i e and g e . In or- 



six quantities Sn e , <p, An, i 
der to have a closed system, we need find additional 
equations for two of these quantities. The ion gy- 
rokinetics will give us this additional information: 
namely, we will determine the ion parallel flow ve- 
locity uiu, which turns out to vanish to lowest order, 
and 5n e /no e , which, because of the quasineutrality, 
is the same as the ion density perturbation Srii/noi. 

As we explained in Sec. Ill BL the ions remain fully 
gyrokinetic under our ordering, i.e., the Bessel func- 
tions in the expression for (x) R . [Eq. ([7])] cannot be 
expanded in small argument because aj ~ k^pi is 
finite [Eq. (fTU)) ]. However, a significant simplifica- 
tion of Eq. for ions is possible because, as in the 
case of the electrons, the 8B\\ term is an order of 
P e smaller than the ip term [Eq. (|2Tj) ] and also be- 
cause, unlike for the electrons, the An is also small: 
indeed, noticing that UthiA'the = {Tnie/mi) 1 ' 2 and 
comparing with Eq. (|28l). we find 



1/2 



«|| ^11 v thi B SBj_ 



c ck± Bq 

Thus, we have for the ions x = or 
(x) Ri ,k = J oK)Vk- 
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T 

~z 



(47) 



(48) 



Because of the way we ordered the electron 
streaming frequency [Eq. (|15p] and the electron col- 
lisions [Eq. the ion streaming frequency 



1/2 



k\\Vthe 



T 
Z 



%u (49) 



and the ion collisions [Eqs. and ([2"5j) [ are small 
and so the corresponding terms in Eq. (|4} for ions 
(s = i) are negligible to lowest order. What remains 
is the following rather simple equation for the ions, 
devoid of kinetic effects except for the ion FLR: 



0, 



(50) 



where we have introduced a new function 

ZeF 0i 



9i = hi 



0/ 



(51) 



Equation (|50|) is homogeneous and has one very 
straightforward solution: 



9i = 0. 



(52) 



Thus, the ion response under our ordering is essen- 
tially electrostatic and the perturbed ion distribu- 
tion function is [see Eq. (JTJ) ] 



Sfi 



ZeF, 



(Mr, 



(53) 



Note that the first term here is a function of the 
gyrocenter variable while the second is a func- 
tion of the position variable r. When computing 
the ion density and flow velocity, we must integrate 
5fi over velocities while keeping r constant, which 
means that terms dependent on R^ must be gyroav- 
eraged at constant r. Just like the gyroaveraging at 
constant R^, this can be expressed in Fourier space 
as a multiplication by the Bessel function Jo(ai). 
Thus, we calculate 



Srij 



= JL f d \ (6fi) t 
noi J 



(i-r ) 



Zeip 



(54) 



where Tq is the operator that is the inverse Fourier 
transform of 

ToM = — [d 3 v {M ai )} 2 F ai = h(a t )e~ a ^ 
n i J 

(55) 

where a, = k 2 L pf/2 and Iq is the modified Bessel 
function. By quasineutrality, Eq. (|54[) also gives us 
the electron density perturbation: 



5n e Z ~ etp 
— U _ 1 o)7f-- 



(56) 



We have recovered a rather popular standard treat- 
ment of the ion FLRi££ Equation (|56|) is often re- 
ferred to as the gyrokinetic Poisson equation. Note 
that Eq. is consistent with the ordering for 

the density perturbation assumed in Sec. Ill Bl [see 
Eq. |Q3J]. 

It will be useful to remember that for k±pi <C 1, 
Tq (a, ) ~ 1 — oti and so 



- (i - r ) ~ kipi 



(57) 
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in the long- wavelength limit — this leads to a signifi- 
cant simplification of Eq. ([56)) . The long- wavelength 
form of the gyrokinetic Poisson equation is useful 
when one wishes to neglect the ion FLR effects. For- 
mally, this can be done by means of the cold-ion 
approximation. r«l. In view of Eq. (|10j) . it imme- 
diately implies k±pi ~ ^fr <C 1. Note that taking 
this limit does not require us to order r within the 
/? e expansion as r and /3 e do not interfere with each 
other in most of the ordering arguments of Sec. Ill Bl 
The only exception to this is Eq. ([24]) . which im- 
plies that if r is excessively small, ion collisions (and, 
therefore, ion viscosity) might become important. 
The condition for this not to happen is f3 e <C t 3 , 
or r » (me/mi) 1 / 3 . Clearly, it is safest to keep 
t only moderately small and to treat any possible 
expansion in r as subsidiary to the /3 e expansion. 

Finally, since 5fi given by Eq. (|53[) is even in vu 
and gyroaveraging only involves v±_ , we have 



1 

not 



d 3 v«|| (6fi] 



0. 



(58) 



This means that under the ordering we have adopted 
the parallel current is carried predominantly by the 
electrons [see Eq. (|57|) ]. 

Equations (JSHJ) and (JSHJ) together with Eqs. (|3"9T 
42]) constitute a complete set. We will assemble and 
summarize them momentarily, but first let us work 
out the collisional terms in our equations. 



D. Collisions 

In order to make our equations useable, we must 
calculate the collision terms. We would like to do 
this by employing a maximally simplified collision 
operator rather than the quantitatively correct but 
cumbersome Landau one. In order to choose the 
most appropriate model, we notice first that the elec- 
tron kinetic equation does not involve any nontriv- 
ial evolution of the v± structure in the distribution 
function. We therefore believe it is acceptable to 
treat the kinetic equation as one-dimensional in ve- 
locity space — formally, this can be thought of as 
integrating out the v± dependence in g e and -Fbe- 
The simplest one-dimensional collision operator is 
the Lenard-Bernstein onej^l which we write in the 
following form: 

dhe 
dt 

2dV\\ \dV\\ V Vthe 



+ ( 1 - 2 ^) F -i/^( 1 - 2 ^ 2 )^ 



, (59) 



where uii = Wii/t^he an d the additional, non- 
differential terms in the operator have been con- 
structed in such a way that electron momentum 
evolution is captured correctly, particle number and 
parallel kinetic energy are conserved and the H the- 
orem is satisfied (cf. Ref. HH; the H theorem will 
be proved explicitly in Sec. IV B|) . The collision fre- 
quency v e i is taken to be velocity-independent — a 
gross simplification compared to the uncensored re- 
ality, in which it is in fact a strong function of veloc- 
ity. We do not believe this to be an important short- 
coming of our model, however, because our focus will 
be on the role of the collision operator as small-scale 
regularization in velocity space rather than on quan- 
titatively precise calculations of the effects of finite 
collisionality. 

Let us now calculate the collisional terms in 
Eqs. QD) and (g2J}. Using Eq. ([33) for h e and 
Eq. ([37| for u\\ e , we get 



1 



dhe 
dt 



dlV\A h (60) 



which simply gives rise to a resistive diffusion term 
in Eq. ([4"0") with Ohmic diffusivity r\ = v e id 2 e - The 
collision operator defined by Eq. (|4"3|) becomes 



C[9e 



1 d 

2di^ 



_d_ 

dii« 



2v, 



1 - 2v\ 



Ik 



Toe 



?0e 



(61) 



where we have used Eqs. ([33]) and (j60|) . the proper- 
ties of g e [Eq. (|35|) ] and the definition of the parallel 
electron temperature perturbation [Eq. (|41|) [. 



E. Summary of the equations 



Assembling now Eqs. ([3JH42]), ([56 [) . (|58 |) and (j60| . 

we arrive at the following closed set 



(62) 



d Z , 

dtT^ 1 


f ) elp - 

- l Q)7f- — 


b- V— d 2 e y 2 { A« 
cm e 11 


d (A 
dt ^ ^ 


-d 2 e VlA {l ) 


- »?V ± A|| c ^ 


cT 0e 

e 


b- v 




-LOe -LOe 


dg e . 


l,b-V 


(g e - 


-^F 0e ) =C[g e 




S) 

V the ) 


F 0e b 


■ V— d^VlAu, 
cm e 



(63) 



(64) 
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where the following short-hand notation is used 



1 



3 2v i 

fv g e 
n 0e J < e 



Toe 



d 8 C r i 



(65) 

(66) 
(67) 



V = v e id\ is the Ohmic diffusivity and C[g e ] is 
the collision operator — a simple model for which, 
based on the Lenard-Bernstein operator, is given by 
Eq. ([CTjl . These equations evolve three fields: tp, An 
and g ei which are all functions of time and three 
spatial coordinates; g e is also a function of v || and 
vj_, although the v± dependence can be ignored (or 
integrated out) if the Lenard-Bernstein collision op- 
erator is used. 

Equations (p2"HM)) constitute a minimal physically 
realizable paradigm for magnetic reconnection and, 
more generally, low-frequency nonlinear plasma dy- 
namics with a strong guide field, including all effects 
we expect to be important: ion sound scale physics, 
ion FLR, electron inertia, electron collisions, Ohmic 
resistivity, and electron temperature perturbation 
determined by a kinetic equation. We will refer to 
these equations as Kinetic Reduced Electron Heating 
Model (KREHM). Its hybrid fluid-kinetic nature and 
the presence of the kinetic electron heating channel 
(further discussed below), constitutes the main dif- 
ference with previously considered models: in the 
2D, collisionless case, our equations can be manip- 
ulated into a form similar to that proposed by de 
Blan k 53 i 54 ; by setting g e = (isothermal-electrons 
closure), we recover the equations of Schep et alM- 
While the use of a fluid model (g e = 0) might be a 
useful simplification, we stress that setting g e = 
cannot be rigorously justified (at least in the ana- 
lytical framework we have chosen): indeed, g e = 
is not a solution of Eq. (|64p unless b • VV^_^4|| = 

(i.e., b • Vjii = 0), which cannot be the case in a 
reconnection-relevant solution. 

We will show in Sec. IIVI that in the collisionless 
(or, more precisely, weakly collisional) limit, the 
coupling of the fluid system to the kinetic equa- 
tion (|64|) via parallel electron temperature fluctu- 
ations provides the electron heating channel and so 
makes collisionless reconnection thermodynamically 
irreversible. We believe this to be fundamentally im- 
portant and physically the most interesting outcome 
of our calculation. In the remainder of the paper we 
will concentrate on this aspect. 

We stress that although the kinetic coupling leads 
to irreversibility and heating, it does not constitute 
a flux-unfreezing mechanism by itself; magnetic field 
lines can only be broken by resistivity and electron 



inertia (a simple proof of this statement was pro- 
vided in Sec. MM- 

To conclude this brief summary of our equations, 
let us note that a number of well-known limiting 
cases are easily derivable from them. As already 
discussed, in the collisionless limit, for g e = 0, they 
reduce to the two-fluid model for strong-guide-field 
collisionless reconnection^ In the collisionally dom- 
inated limit, Eq. (fM| reduces to the standard equa- 
tion for the evolution of the parallel electron tem- 
perature via parallel heat conduction (see Sec. lVCl) . 
If collisions are so large that the resistive scale is 
larger than both p s and d e , Eqs. (|62|) and (|63|) reduce 
straightforwardly to the standard Reduced MHD 
equations^ In the intermediate limit k±pi 3> 1, 
k±d e <C 1, they reduce to the low-beta limit of 
the Electron Reduced MHD equations,— and so can 
support kinetic Alfven waves. More generally, we 
show in Appendix [B] how the full collisionless gy- 
rokinetic dispersion relation in the asymptotic limit 
of k±p e <C 1 and low betaS. is recovered from 
Eqs. ([52TIM)) . Finally, for the linear tearing mode, 
the paradigmatic linear problem of magnetic recon- 
nection theory, Eqs. (|62If6~i| recover the correct ki- 
netic formulation both for the collisionless and the 
semicollisional regimes^ — this topic is also treated 
in detail in Appendix [Bl 



IV. ENERGETICS 

A. Free energy 

A broad class of Sf kinetic systems, including 
gyrokinetics, conserve (in the absence of collisions) 
a positive-definite quadratic quantity that has the 
physical meaning of the free energy of the com- 
bined system of particles and perturbed fields and 
plays the role of the generalized energy invariant (see 
Refs. |48||49| and references therein): 



W 



E 



d 3 r 
~V 



2F, 



Os 



d 3 r \SB\ 



V 8tt 



= -J2 T 0sSS s + U, 



(68) 



where SS S is the perturbed entropy of species s and 
U is the energy of the perturbed magnetic field. 
Since energy flows play a fundamental role in all 
nonlinear phenomena, it will be instructive to under- 
stand the energetics of our equations. As KREHM 
is a partucular limit of gyrokinetics, this is done by 
specializing from the more general gyrokinetic case. 
We saw in Sec. IIII Al that 

*/.= f^ + ^W + * (69) 



n-Oe 



"the 



11 




[see Eq. ([3"3"]) [. Hence we find immediately that the 
electron perturbed entropy is 



— To e SS e — 



The three terms here are the electron density vari- 
ance (denoted Y and, as we shall see, interpretable 
as generalized enstrophy) , the kinetic energy of the 
parallel electron flow (denoted K e ) and the free en- 
ergy associated with the reduced electron distribu- 
tion function <? e , which we will refer to as electron 
free energy and denote H e . Using Eqs. ([37]) . ([56)) 
and ([55} , we rewrite the above expression as follows 



T 0e SS e =Y + K e 



k 



[1 - r (aO] 



2 e 2 |y k | 2 



d 3 rd 2 \V 2 An 



V 



87T 



d 3 r 



3 . ^Oeffe 



2F, 



Oc 



(71) 



where To was defined in Eq. (|55|) . 

Using Eq. ([53"]) for <5/i, we find the ion perturbed 
entropy: 



ToiSSi 



- [i -r (ct;)J — 

1 0e 



(72) 

Note that, were this quantity restricted to kxpi <C 1, 
it would simply be the kinetic energy of the E x B 
flows: indeed, using Eq. (j57[) and assuming an ex- 
pansion in t (as explained just after the latter for- 
mula), the quantity given by Eq. (|72p becomes 



K, 



d 3 re 2 n 0e 2 2 
-Tr^p r —Ps\V±<p\ = 



d 3 r 
~V 



(73) 

where u± = c\V±(f\/Bo is the E x B flow velocity. 
Under the same approximation, the first term on 
the right-hand side of Eq. ([7T|) (which arose from 
the electron density variance) is recognizable as the 
enstrophy of the E x B flow: 



Y 



d 3 re 2 n 0e 4 2 2 



V 2T 0e 
Finally, the magnetic energy is 

d 3 v5B 2 fd 3 r\V ± A«\ 2 



U 



V 8tt 



V 8tt 



(74) 



(75) 



because 8Bu is subdominant under our ordering 
[Eq. pi])]. Combining Eqs. (HE]), G2) and (JTSJ to 



reassemble the total free energy [Eq. ([55]) ]. we get 



W 



K l + Y + U + K e + H, 
Z 



E 



i 



(i-r ; 



Z ft r ^ ^Oel^kl 2 

— I 1 - J-o) — ™ 

T 2J 0e 



d 3 r |V ± A||| 2 + d 2 |ViA||| 2 



y 


8tt 








2-TOe 



(76) 



In three dimensions (3D), W is the only quadratic 
invariant of our equations. In two dimensions 
(2D), there is an additional family of invariants 
(inherited from the more general 2D invariants of 
gyrokinetics^) — they are worked out in Appendix 
|A"1 Their existence, while opening interesting av- 
enues of investigation of an academic kind, suggests 
that one should be very cautious in generalizing 2D 
analytical and numerical results to 3D reality (sim- 
ilarly to the situation in fluid turbulence theory, 
where nonlinear interactions and energy flows are 
dramatically different in 2D and in 3D). 



B. Energy exchange between fields 

It is illuminating to consider separately the time 
evolution of the five constituent parts of W . Multi- 
plying Eq. ([62]) by ep/T 0e or by (Z/r)(l-f )e^/T 0e 
and integrating over space, we obtain the evolution 
of the ion kinetic energy and enstrophy, respectively: 



dK, 
dt 



dY 
dt 



c f d 3 r - „ - A 



d 3 r 



T J ' b " v 



To)<P 



(77) 
(78) 



where Eft = —h-Wip is the electrostatic part of the 
parallel electric field and jy = — (c/47r)VjAi is the 
parallel current. Multiplying Eq. (J63j) by V 2 L A||/47r 
and integrating over space, we obtain the evolution 
of the combined magnetic and electron kinetic ener- 
gies: 



dt 



-J||b-V 



Ke) = 

z 



V 



4tt 



(i-r )v 



i 



j||b-V<JT| 



Ik 



(79) 



With the help of Eqs. (|77H75]) . we now exam- 
ine the evolution of the combined "fluid" (electro- 
magnetic) part of the free energy, Wfl u id = K * + 
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Y + U + K e , which is the quantity that is normally 
considered to be conserved in two-fluid models of 
Hamiltonian collisionless reconnection^! The en- 
ergy exchange terms containing ip cancel and we 
find that Wfl u id can only change due to two effects: 
the resistive Ohmic dissipation [the first term on 
the right-hand side of Eq. (|79|) ]. and an exchange 
with the electron free energy H e controlled by the 
last term on the right-hand side of Eq. ([75)) . Since 
j|| = — en 0e ti|| e [see Eqs. (|37[) and (|58[) [. we can inter- 
pret this term as work done by the parallel electron 
pressure 5p\\ e = no e 5Ti\ e : 

f d 3 r 1 - r f d 3 r 

/ ITe • 7|t l|e = ~ / — U \\eh-V6p lle = -Q. 

(80) 

Thus, the fluid part of the free energy evolves ac- 
cording to 

d „ 47r f d 3 r 9 . . 

5 ^ = -fl-^/ yJ i. (81) 

Note that Q is not sign-definite and can correspond 
both to loss and gain of energy. In particular, it 
clearly represents a loss (Q > 0) if the electron fluid 
is compressed. 

The fluid energy lost or gained this way is re- 
covered in the evolution of the electron free en- 
ergy, which we obtain by multiplying Eq. (|64j) by 
Toe<7e/foe and integrating over the entire phase- 
space (positions and velocities): 

= Q - D colh (82) 

dt 

where -D co ii is the dissipation term due to collisions, 
which must be positive definite for any collision op- 
erator that satisfies Boltzmann's H theorem (see dis- 
cussion in Ref . HH and references therein) : 

Adding Eqs. (JHTJ) and (fS2"]). we find that the total free 
energy W is conserved in the absence of collisions, 
as stated at the beginning of this Section. 



C. Dissipation, electron heating and 
irreversibility of collisionless reconnection 

Consider some initial configuration prone to re- 
connection. Such a configuration will possess a cer- 
tain amount of magnetic energy U, which in the 
process of reconnection will be converted into some 
other form. There can be two fundamentally differ- 
ent types of such energy conversion. 



First, since the conserved quantity is not U but 
the total free energy W = U + K e + K { + Y + H e , 
the magnetic energy can be transferred dynamically 
into K e , Ki, Y or H e without loss of W, or increase 
of entropy — therefore, in principle, reversibly. 

Second, the magnetic energy can be dissipated via 
the resistive term in Eq. l[8Tj) and/or via the colli- 
sional dissipation term D co \\ in Eq. (|82p — in the 
latter case, it first has to be converted into the elec- 
tron free energy via the energy exchange term Q in 
Eqs. (f8Tj) and (|82|) . These processes are irreversible 
because they involve conversion of the fluctuation 
energy into the thermal energy of the bulk electron 
distribution, or electron heating. Indeed, it is not 
hard to sho w 48 ' 61 that 



3 dT 0e fd 3 r f T 0e 5f e fdSf e \ 

where Eq. (|55|) was used to express 5f e and over- 
bar means averaging over dynamical timescales (in 
gyrokinetics, the rate of change of the mean equilib- 
rium quantities is ~ e 2 w, so the transport equations 
are obtained via intermediate time averagin g 61 ' 90 ). 

If we formally forbade any collisions at all and 
set v e i = exactly, the exchanges between different 
constituent parts of W would be the only possible 
energy conversion mechanism. Indeed, Hamiltonian 
theories of collisionless reconnection based on con- 
servative fluid models (g e = 0, Wfl U id = const) find 
that the nonlinear stage of a reconnecting mode cor- 
responds to a transfer of magnetic energy U into 
ordered parallel electron kinetic energy K ei perpen- 
dicular ion kinetic energy Ki, and enstrophy of the 
E x B flow Y . In this situation, the only way energy 
can be lost is via ejection of material — although the 
reconnection process remains technically reversible. 

Is this, however, a good approximation of what 
will in fact occur? We believe that a priori it is 
not because, even if the collisionality of the plasma 
is small, the dissipation terms cannot be neglected 
as large spatial gradients will imply finite currents 
and large velocity-space gradients will imply finite 
values of C[g e ]. The latter effect is especially impor- 
tant. Indeed, if resistivity (Joule heating) is ignored, 
magnetic energy can still be converted into electron 
heat via transfer into H e , formation of small-scale 
structure in v\\ by means of linear phase-mixing (the 
second term on the left-hand side in Eq. (|64|) : see 
further discussion in Sec. [V| and the consequent 
collisional dissipation of the resulting fine-scale elec- 
tron distribution: since the collision operator is a 
diffusion operator in vu [see Eq. (|6T[) ]. -D C oii will 
have a finite value provided structure develops in 
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the velocity-space that satisfies 

Svu 



52 > 



< 



Vthe 



1/2 



(85) 



This phase-mixing channel of electron heating is 
open even in very weakly collisional plasmas. It 
seems a priori clear that its presence should be qual- 
itatively important as it breaks the Hamiltonian na- 
ture of the problem and renders the process of "col- 
lisionless" reconnection irreversible. 

The question of course remains whether the phase 
mixing will be important in any given physical situ- 
ation of interest. Since most such situations are non- 
linear, it is difficult to provide a universal answer to 
this question beyond the general argument that na- 
ture rarely ignores an energy dissipation channel if 
one is available. Here it is opportune to remind the 
reader that there have been several recent numeri- 
cal studies of nonlinear kinetic reconnection that re- 
ported the non-negligibility of the parallel electron 
temperature gradient term in the generalized Ohm's 
law2^ii££2 [see Eq. flS3J}] — and, therefore, the 
non-negligibility of the energy exchange term Q in 
Eq. (fSTj) , which ultimately transfers free energy into 
electron heat [Eq. (|82p]. It is a testable prediction 
that in all such cases, small-scale structure should 
be generated in the velocity space (see Sec. lVD[) . so 
as to enable the electron heating channel. 

In Sec. [VI we discuss the electron kinetics a little 
further, in particular quantifying the notion of the 
free-energy transfer to small parallel scales in ve- 
locity space. We will conclude (in Sec. IV Ep that for 
configurations that evolve sufficiently slowly in time, 
the electron heating rate is finite and independent of 
the collision frequency as v e i — >• +0. 



V. VELOCITY-SPACE DYNAMICS 
A. Hermite expansion 

The velocity-space dynamics and the emergence 
of small-scale structure in v\\ are best understood in 
terms of the expansion of the electron distribution 
function g e in Hermite polynomials i^Ir— Let 



^2, H m (v\\) 



g m F 0e (v\\), (86) 



where = v^/vthe, perpendicular velocity depen- 
dence is understood to have been integrated out, the 
Hermite polynomials are 



H m (v\\) = (-l) m e*ii 



dvl; 



(87) 



(so Ho = 1, Hi = 2u||, H2 = 4t)jj— 2, etc.) and <? m are 
the Hermite expansion coefficients (dimensionless), 
which can be calculated according to 

9n 



dVH 9e(V\\). 



noe J -00 " V2 m m! 

In view of Eq. (|36j) . we must have go = .9i = 
Equation (|4ip implies, by definition, that 

1 ST\\ e 



<n = 



(89) 



y/2 Toe 

Using the orthogonality of Hermite polynomials, 

/ dv \\ ^7— j F 0e = S mn , 90 

n 0e J-oo 2 m m! 

we note that the electron free energy in terms of the 
Hermite coefficients is 



H„ = 



V 



Or 



Toe9 2 e 
2F, 



d 3 r no e To e ^-v ^ 2 

~V 2 ^ 9m 



(91) 



Taking the Hermite transform [Eq. (|88|) [ of 
Eq. (|64|l and using the recursive property of the Her- 
mite polynomials, 

0||.ff m (v||) = - H m+ i(v\\) + mi? TO _i(t)||), (92) 

we arrive at 
dg m 
dt 



Vtheh ■ V 



lm + 1 



9m-l 



cm, 

Vei (mg m - 2<5 rn,2 92) 



d 2 e VlA l{ 



,192 



(93) 



Note that Hermite polynomials are eigenfunctions of 
the Lenard-Bernstein operator ([61]) . 

It is not hard to see that Eq. (|93|) is consistent 
with go = gi = 0. We would like to recast the 
equation for gi in terms of the electron temperature 
perturbation with the aid of Eq. 



d 5T\\, 

dt Toe 



+ v the h ■ Vv^ 
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-2b - V d 2 e W 2 ± Au. 

cm e 



(94) 

This equation is coupled to the rest of the kinetics 
via the heat flux proportional to 173 (the second term 
on the left-hand side). For m > 3, we have 



dg„ 
dt 



+ Wtheb • V 



-v ei mg m . 



Im+l 



2 ' 9m+l + d 9m-l 



(95) 
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Equation ([9"5)l shows that Hermite modes of order m 
are coupled both to lower (m — 1) and higher (m + 
1) modes. Collisions provide a cutoff at sufficiently 
large m regardless of the magnitude of the collision 
frequency. We will discuss this further in Sec. IV Dl 

Together with Eqs. (JB2) and Eqs. (JM]) and 

(f9l>)) can be solved as a system of partial differen- 
tial equations in the three- or two-dimensional posi- 
tion space. This appears to be an attractive way 
to carry out numerical simulations of collisionless 
(in fact, weakly collisional) reconnection or, indeed, 
of other kinetic phenomena in strongly magnetized 
plasmas. This system is substantially simpler than 
the full kinetio2i~— or gyrokinetici^r— descriptions 
that have been used to study reconnection or plasma 
turbulenc o 52 i 96 so far (see Sec. IV Dl for estimates of 
how many Hermite modes must be kept for any given 
collisionality) . The first numerical study using our 
equations is reported in Ref. 1971 

The Hermite formalism allows us to provide very 
concise derivations of three important results: the 
H theorem for our collision operator (Sec. IV Bj) . the 
so-called semicollisional limit of our equations, in 
which collisionality dominates and the parallel elec- 
tron temperature is determined by a fluid equation 
(Sec. IV Cp . the Hermite "spectrum," which quanti- 
fies the fine structure in the velocity space caused 
by the parallel phase mixing (Sec. IVD[) . and finally 
the electron heating rate (Sec. IV E[) . 



B. H theorem 



As we mentioned in Sec. MI D| a key requirement 
for choosing a model collision operator is that it sat- 
isfies Boltzmann's H theorem, i.e., that collisional 
dissipation leads to increase of entropy. This means 
that we must have -D C oii > [see Eq. Substi- 
tuting Eq. (|86|) . using the collision operator given by 
Eq. (|61[) and the orthogonality of Hermite polyno- 
mials [Eq. 1(90])]. we get 



1 



g e C[g e 



n 0e 7-oo " Foe 
oo 

= -Vei m 9m. < 0, 

m— 3 

so Dcoii > 0, q.e.d. 



(96) 



m=0 



C. Semicollisional limit 

Consider the semicollisional limit, v e i 3> u> and 
k\\X m ipe "C 1 ("semicollisional" because the perpen- 
dicular microscale effects associated with p s and pi 



are retained, although the electron inertia term in 
Eq. (|63p must be neglected compared with the re- 
sistive term). From Eq. (f9l>)) . it is clear that in this 
limit, the Hermite coefficients get smaller with to: 



fj: 



k\\Vthe 



fell A 



|| ^mfpe 



mv ei 



< 1. 



(97) 



Therefore, Eq. (|9"5|) allows us to express higher-order 
coefficients in terms of the gradients of the lower- 
order ones and, in particular, §3 in terms of §2, i.e., 
the heat flux in term the temperature gradient: 
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1 Uthe b- v STl[ ' 



2^3 v e 



T, 



Oc 



(98) 



Substituting this into Eq. (|94p . we obtain the stan- 
dard equation for the electron temperature: 



- «|| e b-V [ b- V- 



dt T 0e 



Toe 



-2b- V d 2 e V 2 ± A l{ . 

cm e " 

(99) 

where K\\ e = v 2 he /2is e i is the parallel (Spitzer) ther- 
mal diffusivity. In the semicollisional limit, Eq. 
replaces Eq. (|64|) and completes what has become a 
purely fluid system. 

The energetics of this fluid system are a simple 
particular case of the general situation discussed in 
Sec. IIVI The electron free energy is [see Eq. ([9"T]) ] 



U e = 



d 3 r n ae T Qe „ 2 
V 2 92 



d 3 r n 0e T 0e ST 2 



V 



1 0e 



(100) 

From Eq. (|99j) . we get an evolution equation for H e 
equivalent to the general Eq. 



~df 



d 3 r n 0e T ae 
V 2 



b V- 



6T\ 



T 0e 



(101) 



where Q is the energy exchange given by Eq. (|80|) 
and the dissipation term is the same as D co ii defined 
by Eq. ([55)) — this is checked by using Eq. ([9"6"| and 
restricting the sum just to 173, as given by Eq. (|98j) . 



D. Hermite spectrum 



We now return to Eq. (|95p , linearize it and Fourier 
transform in the parallel direction. If we denote 
<?m(fc||) = (*sgnfc||) m <? m (fc||), we can write the result- 
ing equation in the following form 



dg 



= -v ei mg m , 



m+ 1 _ 

g m +i 



~2 9 m -1 



(102) 



15 



which has the convenient property of supporting real 
solutions*!^ If we now define the Hermite spectrum 
as E m = \g m \ 2 /2 = </f„/2, we find that it evolves 
according to 



dE m 
dt 



m+l/2 



m-1/2 j 



2u ei mE m , (103) 



This cutoff supercedes the collisional cutoff (| 106[) if 
m 7 <C 77i c, i.e., if the mode is growing so fast that 



7»(|fc||kth e ) 2 M 3 - 



(109) 



Otherwise the steady-state spectrum (|105|) is recov- 
ered. 



where T m _ 1/2 = \k\i\v t bey/rn/2 g m 9m-i can be 
thought of as the flux of the electron free energy 
through the Hermite space. When m 1, we may 
approximate r m w |fcii \vtheV^ m E m and 



dE„, 
dt 



\k\\\Vfhe 



dm 



2vpjmEr, 



(104) 



Hermite spectrum in steady state 



Equation (|104p has a steady-state solution-^ 



E n 



exp 



m 



3/2' 



(105) 



E. Electron heating rate 



The results of Sec. IV Dl allow us to make an esti- 
mate of the electron collisionless heating rate. Using 
Eqs. ([53")) and ([9^)1 . we may write 



(110) 



£>coii = -n 0e T 0e ^2 ^ v &i E mEv 



Approximating X)m=3 TO ^m ~ /q°° dmmE m and us- 
ing Eq. (|105p . we find for the steady-state or slowly 
evolving modes, 

£>coll ~ -no e Toe22 V2\k\\\v t heC(k\\) (111) 



where C(/c||) is some function of fen (determined by 
the dynamics of (p, A|| and 6T^ e ) and the collisional 
cutoff is 



3 |ft|||«the 
2^2 Vei 



2/3 



(106) 



When m <C m c , we have the spectrum scaling as 
E m ~ m -1 / 2 in Hermite space. This scaling im- 
plies that the electron free energy H e cx J2m=2 
is dominated by the Hermite modes m ~ m c , which 
are transferring their energy into electron heat via 
collisions (see Sec. IV E|) . 



2. Hermite spectrum for growing modes 

Consider now a slightly more general situation, 
in which the entire Hermite spectrum is growing 
at some rate 27 (this would be the case, for exam- 
ple, for the tearing mode; see Appendix IB 3|) . Then 
dE m /dt = 2jE m and the solution of Eq. (flDi]) is 



E n 



exp 



1/2 



771 



3/2 



(107) 

where another cutoff has appeared, associated with 
the growth rate 7: 



|fciiKhc 
2V27 



(108) 



independently of the collision frequency v e i as long 
as m c 3> 1, i.e.. v e { <C [fcii [ftixe- Thus, the elec- 
tron heating rate is finite in the limit u e i — > +0. 
This is analogous, e.g., to the situation in standard 
hydrodynamic turbulence where the rate of viscous 
dissipation has a finite limit at small viscosities be- 
cause larger spatial gradients emerge^ — or in sub- 
Larmor gyrokinetic turbulence, where a similar pro- 
cess involves large gradients both in space and in 
v± , 48,49,59 rp ne difference j s that here the mechanism 
for the generation of small scales in uii is linear (par- 
allel particle streaming) rather than nonlinear (tur- 
bulent cascade). 

In contrast, for modes growing fast enough to sat- 
isfy Eq. (|109[) . a similar calculation using Eq. (|107[) 
gives 



D r 



E 

fen 



Vei / |fc||Khe 
2V2 V 7 



C{ku 



(112) 

In the limit i/ e , — > +0, -D C oii — > 0. Thus, linearly 
growing modes do not produce any electron heating 
in the collisionless limit, provided the growth rate 
continues to satisfy Eq. (|109|) in this limit. 



VI. DISCUSSION 

In the above, we have already provided the 
motivation behind this work as well as a non- 
mathematical preview (Sec. fl}, a summary of our 
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equations (Sec. IIII E[) and the general arguments 
leading to the conclusion that electron heating must 
be an important dissipation channel (Sec. lIVCj) . We 
developed this last topic somewhat further in SecfVl 
by setting forth what we believe to be a rather use- 
ful practical prescription for further numerical in- 
vestigations, a conceptual view of the phase-space 
dynamics as a cascade in Hermite space, and a sim- 
ple argument implying that the electron heating rate 
remains finite in the limit of vanishing collisionality 
(except for fast-growing linear modes). 

We do not of course claim that our equations con- 
stitute a general theoretical framework for all types 
of magnetic reconnection. The assumption of low 
beta restricts their applicability only to very strongly 
magnetized plasmas (in Sec lIIBTl we indicated that 
typical real-world plasmas where beta is suitably low 
include the solar corona^ the LArge Plasma Device 
at UCLA 8 - 1 - and edge regions in some tokamaks 8 ^). 
We do, however, believe that at least qualitatively, 
our equations might go beyond their formal domain 
of validity and capture much of the essential physics 
of gyrokinetic reconnection. We do not know how to 
derive a similarly simple and yet physically realizable 
model for high-beta and low-guide-field situations, 
which occur frequently in astrophysical contexts. 

An unresolved limitation of our model concerns 
the problem of ultrasmall perpendicular scales. It 
has been found in several (2D) theoretical and nu- 
merical studies of various two-fluid models of colli- 
sionless reconnection that sub-d e structures can form 
in the reconnection region, giving rise to unbounded 
gradients i 58 i 80 i 99 i 100 The only known way of tam- 
ing these singularities is via hyperdiffusive terms in 
the generalized Ohm's law, provided, e.g., by elec- 
tron viscosity. Since we do not have these terms 
in our equations, it is not guaranteed that the sin- 
gularities can be dissipated via the electron heat- 
ing channel* 1 ^ If they are present, any practical 
numerical study of our equations will require addi- 
tional hyperdiffusive regularization — a somewhat 
embarassing solution, since we would ideally like to 
be able to attribute all field-line breaking and elec- 
tron heating to physical mechanisms that are legit- 
imate under our ordering assumptions. We should 
like to note, however, that the evidence of singular- 
ities stems from 2D analyses of Hamiltonian recon- 
nection models. This means that the phase space 
of the system is constrained both by the Hamilto- 
nian structure and by the presence of multitudinous 
additional 2D invariants (see Appendix [A| , so it is 
not, in fact, inevitable that these singularities will 
still present in a 3D kinetic situation. 

It seems clear that the next logical step is a nu- 
merical investigation of the solutions of our model. 
The key physical questions are whether reconnection 



is indeed made irreversible by the electron heating 
channel, how much energy is converted into heat, 
whether the reconnection is fast and what the struc- 
ture of the reconnecting region is (X point? sta- 
ble current sheet? unstable current sheet?). It 
would also be important to find out whether the 
equations are well posed without additional perpen- 
dicular regularization — i.e., if singularities form. 
All of these issues are likely to be tied up with 
the question of whether the 3D case behaves dif- 
ferently from 2D and is, therefore, the right case to 
investigate. Finally, from a technical point of view, 
our model requires fairly modest computational re- 
sources, so it may prove to be a nimble tool for scout- 
ing out the parameter space in preparation for the 
more computationally demanding gyrokinetic recon- 
nection studies i 46 ' 97 
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Appendix A: Two-dimensional invariants 

1. Lagrangian conservation properties in 2D 

Whereas the free energy [Eq. (|68|) ] is the only 
known invariant of gyrokinetics in three dimensions 
(3D), formally restricting consideration to two di- 
mensions (2D) introduces a host of new conservation 
properties. Indeed, consider the gyrokinetic equa- 
tion Eq. ([4| and with d/dz = and no collisions. It 
can then be written as an advection equation of a 
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single scalar quantity: 



^ + ^{(^,^4=0, (Al) 



7 , Is (X) R . 
">s — IT'S TP, fQs- 
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O.s 



(A2) 



This means that the volume integral of any function 
of /i, is conserved: 



d f d 3 R s 



dt. V 



f(h s ) = 0. 



(A3) 



(this is equivalent to picking vu = iuthe/v^ and 
g e = in Eqs. (|A^|) and (|A5jO . Note that 

c\/2/vthe = cd e /vAPs- It is perhaps useful to 
give the evolution equations for A± with three- 
dimensionality, Ohmic resistivity and nonisothermal 
electrons retained, to show how the Lagrangian in- 
variants are broken by all of these effects: 



d_ , Vthe d 
dt V2 dz 



A-. 



B, 



b- V-STu, 

e 



(A8) 



Note that h s is a function of v\\ and v±, so the above 
equation defines an infinite set of invariants both in 
the sense of the arbitrariness of the function / and 
in the sense that the invariants are parametrized by 
the velocity variables. 

Recalling our solution for the ion distribution 
function (Sec. IIII"C|) . we find immediately that hi — 
gi = 0, the conservation law Eq. (|A3[) is satisfied 
trivially. For the electrons, using Eqs. (f5^|) . (|3"3"|) , 
([55]). (J5BJ), and flSEJl, we find 

h e = — (A\\ - d 2 e ViA l{ ) F 0e 

--(l-f )^F 0e +5e. (A4) 

This quantity obeys the advection equation (jAip . 
which becomes 



For a model with no collisions, isothermal electrons 
and cold ions, these equations reduce to those of 
Ref. llOll In the latter reference, its authors point 
out that, while the Lagrangian invariance of A± is 
lost in their 3D system, there is a cross-helicity in- 
variant that continues to be conserved. We see, how- 
ever, that this conservation law does not survive the 
inclusion of nonisothermal electrons: 

d_ fd 3 rAl-A 2 
dt 



V 



2^2c 

Vthe 



An 



■ b-v-tfzO -(i-foV, 

e " / t 



(A9) 



another constraint on the evolution of the system 
that disappears in our treatment. 



dh e c ~ 



(A5) 



2. Quadratic 2D invariants 



where x = <P - {v\\/c)A\\ [see Eq. ([25])]. Thus, h e , 
parametrized by v\\ , is an infinite family of "La- 
grangian invariants," each with its own stream func- 
tion (c/ Bo)x(v\\). This is the "foliation" of the elec- 
tron distribution function discussed in Refs. I55ll56l . 

While (as argued in Refs. l55ll56h this may be remi- 
niscent of the Lagrangian properties of the two-fluid 
models setting g e = in Eq. (|A4[) does not return 
the Lagrangian invariants enjoyed by that model — 
those have to be derived separately and are particu- 
lar to the isothermal-electrons closure (which is not 
a surprise because g e = is not a solution of our 
equations). For the record, in our notation, they 



ari 



,37 



A± = A,, - ^Vi^ll ± — -(l-fote, (A6) 



"the T 

advected by the effective potentials 

Vthe , 

^ = ^^A|, 



(A7) 



It is usually the quadratic invariants that are of 
most practical interest in nonlinear dynamics (or at 
least they are the ones that have the most physi- 
cally transparent consequences, e.g., the inverse cas- 
cades in 2D gyrokinetic turbulence^) . In gyroki- 
netics, specializing to the quadratic function / in 
Eq. (|A3|) . gives one the following family of 2D gy- 
rokinetic invariants ^S, 



Is (x) 



R 



-F 0l 



V 2Fq s \ Tq s 

(A10) 

We have also now integrated over the velocity space 
to produce a conservation law that globally con- 
strains the system in the phase space. Using 
Eq. (|A4[) to work out the electron invariant, we get 



I e = Y+^X + 2U + K e + H e 
di 



= W + —X + U-K, 



(All) 
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where the definitions of, and evolution equations for 
the total free energy W, enstrophy Y, magnetic en- 
ergy U, electron kinetic energy K e , electron free en- 
ergy H e and ion kinetic energy Ki can be found in 
Sec. IIV Al and the quantity 



X 



does not have a commonly agreed name, but is some- 



times referred to as "Ajj-stuff". It is conserved 
standard 2D magnetohydrodynamics (MHD). The 
2D conservation law we are about to derive is a gen- 
eralization of this MHD result. 

Equation (|A11[) tells us that the electron 2D in- 
variant I e contains the total free energy W, which is 
itself conserved both in 2D and in 3D. The remain- 
der of I e is, therefore, a 2D invariant on its own. 
After straightforward algebra analogous to the cal- 
culations in Sec. IIV 51 we find the evolution equation 
for this quantity: 



d_ 

Tit 



X 



c f d 3 r dA\ 

d 2 e W 2 ^ - - ST l{ 
e " 



-rjU 



-(1-fo) 

T 



(A13) 



This is a 2D collisionless invariant because the 
last term on the right-hand side vanishes if (and, 
generally speaking, only if) d/dz = 0. 

We will not delve any further into the mathemati- 
cal consequences of the 2D conservation laws and in- 
stead limit ourselves to remarking that the presence 
of so many constraints on the dynamics particular 
to the exactly 2D case ought to make one beware of 
too much optimism about the relevance of 2D results 
and intuitions to 3D dynamics. 



Appendix B: Linear theory 

1. Linearized equations 

Since we will want to use the linearized equations 
that are about to be derived for treating the tearing 
mode problem, we would like to do the linearization 
around an equilibrium containing both the guide 
held Bo = BqZ and some in-plane field — gyroki- 
netically speaking, this in-plane field is part of the 
small perturbation of the guide field. Thus, we let 



B = B z + 6B^(x)y + 6B { l 



(i) 



(Bl) 



where 5B^ 
A\°\x)+A 



< SB'y 



(0) 



(i) 
II ' 



« B Q . 

(0) 



Therefore, Au = 



where 5 By ' = —dA^ 
z x VAj 1 '. 



(0) 



/dx = B Q f(x) 

and 5B ( l ] = -z x VA^ 1 '. The function f{x) will 
contain all the information about the in-plane equi- 
librium. Examining Eqs. (JS2HHII)- h is n °t hard to 
see that an equilibrium of this form can be main- 



( A12 ) tained by letting ip^ 



0, g 



(0) 



and adding 



an equilibrium parallel electric field to Eq. ([(k?]) to 
maintain the equilibrium magnetic field against the 



m Ohmic resistivity: cE 



(0) 
I 



r)B Q f'(x). 



Inserting all this into Eqs. (1621 [641 . dropping the 
superscripts on the perturbed quantities and Fourier 
transforming with respect to y, z and time, we get 



Z 



-(l-t )^ = —dl [A,|(*)ViA|| 
t l 0e cm e 

-kyf"(x)A\{\ , (B2) 

- LO [4 1| - rfe V i^||] = -^ViA|| + kyCd 2 J"{x)ip 



— k\\ (x)c 



Lug e + fc||(x)w|| 



Z 



(i-r 



<p-- $ T \\. 

e 



(B3) 




STu 



T, 



Qc 



F 0e = -iC[g e ] 



1 ~^ L ) F ^ d2 e [*l|(*)ViA|| 



-k y f"(x)A l{ ] , (B4) 
where ki\(x) = k z + ,f(x)k y and = d 2 — k 2 



a. Collisionless limit 

The simplest approach to linear theory in a kinetic 
plasma is to consider the purely collisionless case. 
Although, as we explained above, this is not really a 
good limit, even in a weakly collisional plasma, the 
linear results obtained in this limit are useful in that 
they will allow us to make contact with the exist- 
ing theories. The collisionless limit is also a useful 
route to certain linear results (like Landau damp- 
ing) that do in fact depend on infinitesimal amount 
of velocity-space dissipation. 

Thus, let us set r) = and C[g e ) = in Eqs. (|B3|) 
and ()B4j) . respectively. If we now solve Eq. (jB4|) for 
g e explicitly and then integrate over velocity space 
according to Eq. (p5)l to obtain the parallel elec- 
tron temperature perturbation, we get, after stan- 
dard alegebra, 



Ik 



z(C(x)) + C(x)Z'(C(x)) 



T 0e |&||(a:)|vthe Z '(C( X )) 

x—d 2 e [fc„ (x)Vj_A|| - k y f"(x)A ll ] , (B5) 
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where ((x) = ui/\k\\(x)\vthe, Z(Q is the plasma dis- 
persion functioiJ^ (not to be confused with Z in 
the ion charge Ze) and Z'(() = -2[1 + (Z((,)]. All 
kinetic effects are wrapped up in the above expres- 
sion for the parallel elecron temperature perturba- 
tion. Using Eq. (|B2j) . we can simplify it somewhat: 



l z(0 + cz^)z f 

-JT\\e = ~* ^0 - {1 - T ° )lP 



1 - 2C 



z'(0 



z 



(1 - f oV (B6) 



This can now be substituted into Eq. (|B3[) . where- 
upon Eqs. (|B2j) and (|B3j) form a closed set. 

Recalling the gyrokinetic Poisson equation (f5G]) . 
Eq. (|B6|) can be interpreted as an equation of state 
with ^-dependent effective adiabatic exponent: 



G 



(G-l) 



n 0e 



z'(0 



(B7) 
(B8) 



This notation will be useful in our treatment of the 
tearing mode. 



b. Semicolhsional limit 

In the semicolhsional limit, we abandon Eq. (IB4I) 
and instead use the linearized Eq. to determine 
the parallel electron temperature perturbation: 



Ul + iK\\ e k 2 (x) 



Ik 



Toe 



— d\ [k ll (x)WlA ll -k y f"(x)A ll ] 
cm e 



Using Eq. (|B2|) . we rewrite this as 



(B9) 



2. Kinetic Alfven waves 

a. Collisionless limit 

Let us make a digression and ascertain that the 
equations we have derived contain the wave dynam- 
ics that they are expected to contain. Formally, 
if we assume k z ^> k y f(x), k z ^> k y d 2 f"(x) and 
k z d 2 /dx 2 3> k y f"(x), we find ourselves in a homo- 
geneous plasma (all terms containing /(x) can be 
neglected and fcj| = k z ). This means that we can 
now also Fourier transform in x and, after a few 
lines of standard algebra, we obtain from Eqs. (|B2|) . 
(|B3[) and (|B6|) the following dispersion relation: 



k\d 2 j2 



zi-r„(*ip?/2) 



i 



(z(()} = ^kU 2 



(B12) 

where k\ = k 2 + k 2 . It is not hard to check that 
this agrees with the gyrokinetic dispersion relation 
at low beta derived in Ref. ID [their Eq. (D17)]. 

Looking for solutions with Q = w/|fc|i|uthe *C 1 
and using Z(Q « i^/n — 2£, we find the well known 
dispersion relation for Alfven waves (both MHD and 
kinetic): 



w = ±k\\v A k_ L p l] 



1 



i-r (fcip?/2) 



The MHD Alfven waves, 

ui = ±k\\VA, 



(B13) 



(B14) 



are recovered for k±pi -C 1 (in which case r 



1 - kj_p 2 /2), the kinetic Alfve n waves, 



u> 




(B15) 



for fc^pi 3> 1 (r ~ 0). Finally, the damping rate is 
found from Eq. (|B12[) as a small perturbation of ui: 



- ST n „ = - 



2uj 



Z 



lu + inu e k 2 (x) r 



(1-ToV. (BIO) 



t.2 2 



(B16) 



This can be substituted into Eq. (|B3j) . where also 
the resistive term is now retained, while the electron 
inertia must be neglected in comparison (because 

U e i > Ul). 

Recasting Eq. (|B10[> in the form (|B7]). we find the 
effective adiabatic exponent in the semicolhsional 
case: 



G = l 



2uj 



Ul + iK|| e fcn (x) 



(BH) 



(note that this is valid for any k±pi). This again is a 
well known formula for the electron Landau damping 
of Alfven waves at low beta [cf. Eq. (63) of Ref. [6l| 
noting that /3 e = (Z/r)/3j]. 

Another familiar and useful limit is the long- 
wavelength, cold- ion approximation, k±pi ~ yfr <C 
1 (already discussed in Sec. IIIICp . Using Eq. ((57)) 
in Eq. ((BT2]) . we get 



2^ 



1 



[i + CZ(C)] = 



(B17) 
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Assuming again ( < 1, we get for the Alfven waves 
in this limit: 



to read 



Z 



1 = 



(B18) 



k y f (x)v t he T 



(i-r V = 



k\i\v A ^-\n — — . B19 
2 V m l p e 



-U 

2 e 
k y f(x)v th , 



3 2 .2 /"(*) 



/(*) 



"the 



A,,, (B21) 



The dispersive kinetic Alfven waves are be- 
lieved to play a key role in enabling fast kinetic 



i + G-(i-r ) 

T 



(B22) 
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so it is important that they are where G(£) is the effective adiabatic exponent given 



reconnection 
fully retained in our equations. 

Note that assuming £ <C 1 in the above calcula- 
tion meant that we effectively adopted an isothermal 
limit [see Eq. (|B6|) ] . Enforcing this requirement on 
the frequency [Eq. (|B13[) [ also eliminates the elec- 
tron scales: k±d e <C 1 (in the cold-ion limit, this is 
modified to d e <C p s )- The non- isothermal effects 
will come in at the electron inertial scale together 
with the breaking of the flux conservation. 



b. Semicolhsional limit 

Finally, in the semicolhsional limit, using 
Eq. (|B10|) instead of Eq. (|B6|) , we get the following 
dispersion relation: 



by Eq. ()B8|) . Q(x) = /\k y f(x)\v t he, and we have 
neglected d? e k^ and d\f' '(x)/ ' f(x) compared to unity 
in Eq. ([B22| . 

As in all tearing-mode calculations, these equa- 
tions will turn out to have a boundary layer around 
a: = 0. The reconnection of the magnetic flux will 
happen inside this layer, while outside the dynamics 
will be essentially MHD. The difference between the 
standard MHD situation^ and what we have here 
is that there are two microscales in the problem: p s 
(or pi) and d e , so the boundary layer will in fact be 
made of two nested boundary layers. We will first 
deal with the outer solution, then give a set of a pri- 
ori estimates of the relevant scales and of the growth 
rate of the mode that follow from the structure of 
the problem in the inner region. 



a. Outer (MHD) region 



uj(uj + iryfcjj = k^v A k\p\ 
Z3w + iK\\ e ku 



1 



T Ul + iK\\ e k?, 



i - r (fci P f /2) 
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The isothermal limit, Eq. (|B13|) . this time with resis- 
tive damping, is recovered for K|| e ^| 3> w , or j equiv- 
alent!^ k±d e <C fc||A m f pe . 



3. Collisionless tearing mode 

Let us now consider the limit opposite to that 
which produced a homogeneous situation, namely, 
k z <C k y f(x), so the problem becomes effectively 
2D and fcy(x) « k y f(x). This is the tearing mode 
problem. Since we are expecting a (purely) grow- 
ing mode, let u = i"f. In the collisionless limit, 
Eqs. (|B2[) . (|B3[) and (|B6[) can now be rearranged 



Here we seek the solution at the equilibrium scales 
I, assuming I ~ [f'(x)/f(x)]~ 1 ~ fc" 1 > p s > d e . 
Accordingly, we can neglect d^d^ <C f and (Z/t)(1 — 
f o) ~ - P 2 s {d1—ky) <S 1. We are not assuming that ( 
is either small or large. Under these approximations, 
Eq. (|B22[) becomes simply 



Vthe 



An = ip. 



kyf (x)v t h 

Substituting this into Eq. (|B21|) . we get 

/"(*) 



(B23) 
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A|| = ^CW(^-^)C(x)A||. 



(B24) 

Now note that, since y2p s /d e = Vthe/fM, we have 



V2p s 



V) 



k y f(x)\vA U A y' 



(B25) 



where ujA y = \k y f(x)vA\ = l&yiMyl is the in-plane 
Alfven frequency. Assuming 7 <C ujAy (tearing mode 
is sub-Alfvenically slow), we may then neglect the 
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right-hand side of Eq. (|B24[) and so end up with 



dlA 



7-2 

N 



f"(x) 



A 



(B26) 



This is the standard MHD tearing mode outer-region 
equation. Its solution depends on the particular 
choice of f(x), but we will not need its detailed form. 
Since SB y = — d x A\\ must reverse direction at x = 0, 
the sign of d x A\\ will be different on either side of 
x = 0, so the derivative will have a discontinuity. 
This is quantified by 



A' 



1 



[d x A 



-o ■ 



(B27) 



the tearing-mode instability parameter. It is a mea- 
sure of the singularity developed in the current by 
the ideal-MHD outer solution as it approaches the 
reconnection layer. The singularity is resolved by 
corrections to the ideal Ohm's law (electron iner- 
tia in the case we are considering). The solution in 
the boundary layer around x = will have to be 
matched to the value of A'. 



by Cowley, Kulsrud and Hahm 6 - 7 - directly from the 
Vlasov-Maxwell kinetics: 



Ai 



ip cr 



1 X 
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= -(l-f )& 
r 

d\d 2 x A\\ = -(1-fo)^, 



(B32) 
(B33) 



where the "scaled effective conductivity" is 

1 IS* 



-(!) 



2x 2 



Z' 



(B34) 

Equations (jB32l) and (|B33|) are Eqs. (26) and (25), 
respectively, of Ref. |67], specialized to the case of no 
equilibrium density gradienfeii^ [note that in Ref. WA , 
5 is defined in terms of u, not 7, so to get exactly 
their equations, one must replace 6 — > —16, a — > —a 
and restore (p = — {c/ vt\ le )iip\. These equations de- 
termine the parallel electric field [Eq. (|B32[) ] and the 
parallel current [Eq. (|B33[) [ in terms of the ion re- 
sponse. This response is, in general, nonlocal, which 
leads to considerable technical difficulties ! 66 ! 6 ^ 73 ' 105 



c. Mathematical structure of the problem 



b. Inner region: equations 

At scales of order p s or smaller, we may expand 
the in-plane equilibrium profile: f(x) « xf'(0) = 
x/L s . Then C( x ) ~ i6/\x\, where 



7 



1 



kyV A V2p s 



(B28) 



Since d 2 3> fc^, f "(x)/f(x) at these scales, 
Eqs. (|B2T|) and ([B22|) become 



™(l-f o )0 = i<#&4||, 



1 



dldl) A,, = 



l + G-(l-f ) 

T 



<P, 



(B29) 
(B30) 



where f = f (p|9|/2), (p = -(c/v tile )itp, and the 
effective adiabatic exponent G is given by Eq. (|B8j) , 
which we now rewrite as 



= -2 



1 



5 1 

2 + Z'(i6/\x\ 



x 



(B31) 



This function has limiting values G(0) = 3, G(oo) = 
1 (electrons are isothermal for x ^> 6) and a smooth 
monotonic transition in between. 

In order to benchmark our theory against existing 
treatments, let us remark that Eqs. (|B29[) and (IB 301) 
can easily be manipulated into a form that coincides 
with the collisionless tearing mode equations derived 



Let us introduce a rescaling of the spatial vari- 
able £ = x/6- m , where 8 ln reflects the width of the 
inner solution and will be determined later. Equa- 
tions (|B30[) and (|B29|) can be rewritten as follows 



266i n Z 
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dl 



(1 - f 0)^ = £A|f, 



(B35) 



(B36) 



where primes denote derivatives with respect to £ 
and we have used Eq. ([B36P to express the (Z/t)(1 — 
t Q )<p term in Eq. (|B30[) in terms of A',!. 

Let X (£) = £Aj, - A [{ = e{A\\/0'- Then A\{ = 
x'/£. Differentiating once Eq. (|B35[) and expressing 
Ay in terms of \) we g e t 



_ I _ 1 "in p \ / _ "in A "in -/ 

dt\e + 2p ) x die <py- 



(B37) 



We now need a closed expression for Cp in terms 
of An, which is where the nonlocality of the ion 
response — the operator on the left-hand side of 
Eq. (|B36[) — presents a technical challenge. The 
simplest choice is to avoid it by considering the limit 
of cold ions. It turns out, however, that with no ad- 
ditional trouble, it is possible to accommodate (non- 
rigorously) the hot-ion limit by using the Pade ap- 
proximant for the ion response^ 



Z 



(i-r ) 



l-{l/2)p 2 d 2 x 



(B38) 
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Then Eq. (|B36)) becomes 



Ultralow-beta (one-fluid) limit 



<P e 6 in 



S -II 



§TX"- (B39) 



By integrating this equation once, we work out the 
(p 1 term in Eq. (fB37) : 



X-Xo 



Z5l 



(B40) 



where xo is a constant of integration, we have used 
Pi/2p^ = t/Z, and arranged for the overall coeffi- 
cient in the left-hand side to be unity by choosing 
the inner rescaling width 



6 m = (V2p s 5 



1/2 



(B41) 



It remains to substitute Eq. (|B40[) into Eq. (|B37[) 
and arrange terms neatly. The result is 
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(e 2 + A 2 )x = A 2 , 
(B42) 



(B43) 
(B44) 



and G(x/S) is given by Eq. (fB3Tj) , with x/ 5 = V2ct£,. 

The above equation is subject to matching with 
the outer solution. First, using the definition of \ 
and the property A'^ = x'/£; we can write 



A,, 



-Xo 



(B45) 



This is matched to the outer (MHD) solution (see 
Appendix IB 3 aj) by identifying the limit £ — > oo for 
the inner solution with x — > for the outer solution: 
this gives x(£ -> °o) = 0, xo = ~ %(| 7 ^ °°) = 
— A\\(x — > 0) and, from the condition (|B27[) . 



(B46) 



Equations (|B42j) and ((B46j) constitute an eigen- 
value problem, where A 2 is the rescaled eigenvalue 
and a is the parameter that measures the width of 
the ion layer. Due to the double- layered structure 
of the inner region, the exact solution requires some 
work^&SiiflS ( see Appendix EU, but all the rel- 
evant scalings can be derived in a simple heuristic 
way. While these scalings are known, it is useful to 
summarize them in a unified exposition. 



If a < 1 (p s <C 6), the ion effects are negligi- 
ble and our equations turn into the standard tearing 
mode equations^ 2 - extended to arbitrary A /63 i 64 and 
with resistivity replaced by electron inertiai 71 i 72 

It is easy to derive the scaling of the growth rate 
without an exact solution of the equations. From 
Eq. (|B42| . A 2 - x7£- The width of the X (0 func- 
tion is order unity (i.e., x ~ <5j n in dimensional 
terms), so the integral in Eq. (|B46|I is over an in- 
terval of order unity. Therefore 



A 2 ~ A'Su 



(B47) 



Using the definitions (|B4I|) and (|B43j) . one gets the 
scaling a 66 ! 71 



dtA 1 ' 
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d 2 e A\ 7 ~ kyV A 
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(B48) 

This result breaks down when A' is so large that 
A'Ji n ~ 1 or larger, in which case the magnitude of 
the current is limited by the width of the reconnect- 
ing region, 6- m . The correct scalings are obtained by 
replacing A' by l/5 in in Eq. (|B47|1 : 



A 



1 



(B49) 



(in other words, instead of estimating the current 
as 9 2 A|| ~ (A'/5 in )A\i, one takes d 2 A|| ~ A^/SfJ. 
Hence one gets the scaling a 65 ! 72 



Ps 



d e , 'f ~ kyV A -^. 



(B50) 



As we explained above, all of this is valid provided 
p s <C 8, which amounts to 



dl 



p e ^ « {d e A'f 



for the finite-A' scalings (|B48|) and 



5§ ~ (3 e -L « 1 
di m P 



(B51) 



(B52) 



for the infinite- A' limit (|B50j) . Either condition can 
only be satisfied if /3 e is very low indeed — not just 
of order, but much smaller than m e /mi [obviously, 
this is treated as a subsidiary limit within our formal 
low-beta ordering (|12pj. 



e. Two-fluid case 

If a ^> 1, or p s 3> 8, the problem contains both ion 
and electron scales and is no longer mathematically 
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equivalent to the resistive tearing mode [the terms 
involving a in Eq. (|B42j) cannot be dropped]. The 
width of the integration region in Eq. (jB4Gl) is now 
£ ~ 1/a (or, dimensionally, x ~ S). Therefore, 

A 2 ~ A'<5 in a. (B53) 

This gives immediately the following scalings^r— 

(B54) 

These scalings are correct provided A'5; n < 1, or 
A'p^di^ 3 < 1. At larger A', again, like in Ap- 
pendix |B331 <5in limits the magnitude of the current 
and we replace A' by l/8- m in Eq. (|B53[) : 



A ~ a, 



(B55) 



whence follow the scaling o 73 ' 74 ' 76 



,1/3 2/3 



(B56) 

Checking now the condition a > 1, we hnd that 
the fmite-A' scalings (|B54[) are valid provided 

§ ~ ft ?± » {d e A'f (B57) 

and the "infinite-A"' ones (|B56| hold for 



§ ~ ft ^ » 1. (B58) 



Asymptotic solution for the two-fluid 
tearing mode 



In Refs. IHzl the scalings (|B56)) were derived 
in a more mathematical way via a rather involved 
double-layer matching procedure in wavenumbcr 
space under the assumption of isothermal electrons 
(G = 1). Th ere a re several other derivations (no- 
tably Refs. I73fl06h. also analytically quite cumber- 
some. Here we give a calculation that recovers the 
essential result without assuming isothermal elec- 
trons and at the minimal analytical cost. 

We start by noting that, using Eqs. (|B45[) and 
(IB46[) . we may write A\\ in a form that is automati- 
cally matched to the outer solution: 



A i 



Ao 



i r°° 

1 + x + - A'5 in £ + £ J dz 



(B59) 

Our strategy will be first to solve Eq. (|B42[) in the ion 
region and use the solution Xi m the above formula 



for Au ; then to solve in the electron region and use 
the solution Xe in Eq. (|B45[) ; finally to take the large- 
argument asymptotic of the electron solution An e 
and match it with the small-argument asymptotic 
of the ion solution Aiij. 



a. Ion region 

The ion region is x ~ p s , or £ ~ a ^> 1. In 
this limit, the l/£ 2 and A 2 x terms in Eq. (|B42[) are 
negligible (the latter because we are anticipating the 
ordering A 2 ~ a) and G ~ 1. Introducing the new 
variable z = £/ayl + t/Z = x/p T , where p T = 
p s y/1 + t/Z, we can now write the rest of Eq. (|B42|) 
as follows 







Xi n ) 

Z l 



where we have denoted 



A 2 = 



A 2 



a 2 (l + r/Z) 



< 1. 
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(B61) 



Equation (|B60[) can be solved exactly, subject to the 
boundary condition Xi( z °°) 0: 



Ci - y 



dre 2r 



ds 



(B62) 



where Ci is a constant of integration and the pa- 
rameter zq in the particular integral can be chosen 
at will, with the difference absorbed into the homo- 
geneous part of the solution (the constant Ci). 

Substituting the solution (|B62I) into Eq. (|B59|) . 
taking the limit z — > and keeping only the leading- 
order contributions (in 1/a) in all terms of the ex- 
pansion, we get 



-II' 



-Xo 



-A 2 ln£ 



ay/1 + t/Z 



, (B63) 



where we formally ordered A'5- m ~ 1 and anticipated 
Ci ~ 1 (all of this will be checked a posteriori) . 



b. Electron region 

The electron region is x ~ 5, or £ ~ 1/a <C 1. 
In this limit, the £ 2 % term in Eq. (|B42p is negligible 
and we are left with a homogeneous equation for 
X + 1 = x/Xo- Introducing the new variable y = 
we get 



y 



dy 



Xe = 



2a 2 



(B64) 
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The right-hand side is, in fact, small, and so we 
can solve this equation perturbatively, in powers of 
X 2 /a 2 ~ l/a (cf. Ref.LZJ). We will see that we need 

to do this to second order: % e = x^ + Xe^ + Xe^- 
Our boundary condition will be that the current, 
proportional to x'/y, must be even as y — > 0. We 
can also let Xe^O) = Xc?\ty — without loss of 
generality, so Xe (0) = xi o) (0). 

Integrating Eq. (|B64jl to zeroth order (right-hand 
side = 0), we get 



(o) 

Xe' _ n 

Xo 



(B65) 



where C e is a constant of integration. To first order, 
we obtain, again by direct integration, 



XO 
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A 2 f y duu 

2a 2 e J 1 + (u 2 /2) [G(u) + t/Z] 



Finally, to second order, 
A 4 



(B66) 



(2) 
Xe 



du u 



Xo 4«4 °J Q 1 + (u 2 /2)[G(u)+t/Z] 
dv f v dw w 



o v 2 J l + (wy2)[G(w)+r/ZY 

(B67) 



Substituting the solution worked out above into 
Eq. (|B45[) , taking the limit y — > oo and again throw- 
ing out all terms subdominant in l/a, we find 



'XoCe 



1 + I G A£ - A 2 In £ + I G AA 2 £ In £ 



where A = A 2 /a^/l + t/Z and 



Ia = 



dy^/1 + t/Z /V2 
o l + (y 2 /2)[G{y)+T/Z] 



(B69) 



is a number of order unity. Had electrons been 
isothermal (G = 1), it would have been Iq = ir/2. 



c. The dispersion relation 



We now match the asymptotics (|B68j) and (|B63|) 
term by term. The ln£ term can be ignored as long 
as we carry out the matching for £ 3> 1 — we can 
do this because the electron solution is in fact valid 
beyond the electron region (£ ~ l/c*0 all the way 
through the intermediate region £ ~ 1, or x ~ Si n 
(because we have solved to second order in l/a). 

Matching the remaining three terms allows us to 
determine the two constants of integration and find 



the dispersion relation 
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This indeed has the two limits that we intuited in 
Appendix IB 3 el 



A' Sin < 1 => A = 



A' 5^ > 1 => A = 



A 2 



= ^~ A'Sir 
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(B71) 
, (B72) 



so we now know not just the scalings but also the nu- 
merical prefactors (which depend on the functional 
form of the effective adiabatic exponent G of the 
electrons via the constant Iq)- Equation (|B70|) gives 
a smooth connection between the two limits. 



5. Semicollisional tearing mode 

a. Equations 

To treat the tearing mode in the semicollisional 
limit (22, we must use Eqs. (|B2j), (JB3J) and ([BIO]) , 
again taking kii(x) f=s k y f(x). The outer region is 
still MHD as described in Appendix IB 3 al In the 
inner region, instead of Eqs. (|B29|) and (|B30|) . we 
get (the first equation is unchanged) 
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(B74) 

(B75) 



These equations turn out to have a mathemati- 
cal structure that can be exactly mapped onto the 
collisionless casej^l Let us introduce a new scale 



5 " = d e 
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V2p s \kyVA 
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{k y L s S)- 1/2 , 
(B76) 

where S = VA.L s jr\ is the Lundquist number. Then, 
denoting $ = (S v /S)<p, we can recast Eqs. (|B731^B75|) 
in the following form 
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(B77) 
$, (B78) 
(B79) 
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where a = 2r/K\\ e /d 2 v 2 he . For the collision operator 
we chose in Sec. MI D| 77 = v e id 2 , /c|| e = v 2 he /2i' e i, 
and so a = 1; the isothermal closure would formally 
correspond to a = 00. Comparing Eqs. (|B77hlB79| 
with Eqs. (|B29I [B3Tj) . we see that all results obtained 
for the collisionless case can be converted into analo- 
gous results for the semiciollisional case by mapping 



V 

1 



(B80) 



and using Eq. (|B79|) instead of Eq. (|B31|) for the 
effective adiabatic exponent of the electrons. 

Again, these equations can be manipulated into a 
form derived in Ref. H3- Proceeding analogously to 
the way we did in Appendix IB 3 b[ we get 



.r 

5„ 



-4 1, 



x ^\ 1 x 
— $ cr 



-(l-f )$, (B81) 



-(l-f )$, (B82) 

T 



where the "scaled effective conductivity" is 
l + (a/2)(a;/<g 2 



2 + (3 + a){x/5 v ) 2 + (a/2)(x/Sr,) 4 



(B83) 

Equations (|B8T|) and (|B82|) are Eqs. (76) and (78), 
respectively, of Ref. IdTI again in the special case of no 
equilibrium density gradient (to recover their equa- 
tions exactly, replace 5 — > —iS, 5 n — > 5 n /V2 and 

<7->o-/2). 

Using the mapping (|B80|) . we conclude that the 
general tearing mode equations (|B42[) and (|B46[) . de- 
rived in Appendix IB 3 cl now hold with 

1/2 / \ ^/^ 

6 in = (V2 Ps 5 v ) = L s ( (kyLsS)- 1 ' 4 , 

(B84) 



Vll \k y v A 



3/2 



(k y L s S) 1/2 , (B85) 



a 



V2 5 n V2<5 r; L s \kyV A 



-1/4 



{k y L s Sf' A , 
(B86) 



and G{V2a£) given by Eq. (|B79|) . All the same 
mathematical considerations apply, with the (qual- 
itatively inconsequential) exception that the func- 
tional form of the effective adiabatic exponent is dif- 
ferent. 

We conclude by summarizing the heuristically ob- 
tainable scalings for the semicollisional case — these 
can be read off from the results of Appendices IB 31 
and El with the aid of the definitions (|B84HB86p , 
From this point on, all lengths are normalized by L s . 



b. Resistive MHD (one-fluid) limit 

Proceeding analogously to the case considered in 
Appendix IB 3 dl we assume a -C 1, or p s <C S v , and 
recover the classic resistive-MHD tearing modej£2~— 
First, at finite A', Eq. (|B47|) with the definitions 
(|B84| and (|B85|l gives the scaling^ 

S v ~ A*' s pj\kyS)-*' B , (B87) 
S in ~ A a / 5 (k v S)- 2 / 5 , (B88) 
A'^ 5 {kyS)-^ 5 . (B89) 



kyVA 



The reconnection layer's width is <5j n -C S v (because 
a <C 1) and the condition for the above scalings to 
apply is A'S m <C 1, which translates into 



A' « (k y S) 1/3 . 



(B90) 



When this is broken, one gets the "infinite- A"' scal- 
ing (|B49p , whence follow the scalings^^ 



6,, ~ p-^kyS)-^ 3 , 

Sin ~ {kyS) / , 

2- ~ (kyS)-^. 



kyVA 



(B91) 
(B92) 
(B93) 



The resistive MHD results are valid provided a <C 
1, which imposes an upper bound on the ion scale: 

p s <C A n / 5 (k y S)~ 2 / 5 (B94) 

for the finite-A' scalings and 

Ps « (kyS)- 1 / 3 (B95) 

in the "infinite-A"' limit. 

c. Two-fluid case 

When the ion scale is sufficiently large (a 
1), two- fluid effects become important, similarly 
to the case considered in Appendix IB 3 el Using 
Eq. (|B53j) and the definitions (jB84HB86|) . one gets 
the scalings£&££ 



A'V3^/3 M 

(A' Ps ) 1/6 (fc^)- 1/3 

(A'psf/HkyS)- 1 ^ 



kyVA 



2/3 



These scalings hold provided A'8i n <§C 1, or 

a' « P :"\kySf\ 



(B96) 
(B97) 
(B98) 



(B99) 



At larger A', Eq. (|B55[) must be used, whence follow 
the scalings£^>ZO 

8 V ~ p^l 7 {k y Sr A l\ (B100) 

s in ~ P y 7 (fc^)- 2 / 7 , (bioi) 

~ pTikyS)- 1 ' 7 . (B102) 
In both cases, the width of the reconnection layer 



is <5, ; <C ^in (because a 3> 1), which holds if 

Ps > A' 1/5 (fc y 5)~ 2/5 (B103) 
for the finite-A' scalings and 

Ps » (k y S)-^ 3 , (B104) 
in the "infinite-A"' limit. 



* Electronic address: |a.zoccol @physics .ox.ac.uk| 

1 Electronic address: a.schekochihinl@physics.ox.ac.uk 

1 D. Biskamp, Magnetic Reconnection in Plasmas 
(Cambridge Univ Press, Cambridge, 2000). 

2 M. Yamada, R. Kulsrud, and H. Ji, Rev. Mod. Phys. 
82, 603 (2010). 

3 R. J. Hastie, Astrophys. Space Sci. 256, 177 (1998). 

4 R. Giovannelli, Mon. Not. R. Astron. Soc. 107, 2163 
(1947). 

5 P. A. Sweet, Annu. Rev. Astron. Astrophys. 7, 149 
(1969). 

6 J. Dungey, Phys. Rev. Lett. 6, 47 (1961). 

7 C. J. Xiao, X. G. Wang, Z. Y. Pu, H. Zhao, J. X. 
Wang, Z. W. Ma, S. Y. Fu, M. G. Kivelson, Z. X. 
Liu, Q. G. Zong, et al., Nature Phys. 2, 478 (2006). 

8 J. P. Eastwood, T. Phan, F. S. Mozer, M. A. Shay, 
M. Fujimoto, A. Retino, M. Hesse, A. Balogh, E. A. 
Lucek, and I. Dandouras, J. Geophys. Res. 112, 6235 
(2007). 

9 A. Retino, D. Sundkvist, A. Vaivads, F. Mozer, 
M. Andre, and C. J. Owen, Nature Phys. 3, 236 
(2007). 

10 S. Servidio, W. H. Matthaeus, M. A. Shay, P. A. Cas- 
sak, and P. Dmitruk, Phys. Rev. Lett. 102, 115003 
(2009). 

11 J. F. Drake, M. Swisdak, H. Che, and M. A. Shay, 
Nature 443, 553 (2006). 

12 J. F. Drake, M. Opher, M. Swisdak, and J. N. 
Chamoun, Astrophys. J. 709, 963 (2010), 0911.3098. 

13 D. A. Uzdensky, Space Sci. Rev., in press, e-print 
larXiv:1101. 24721 

14 A. Bhattacharjee, Annu. Rev. Astron. Astrophys. 
42, 365 (2004). 

15 P. A. Cassak, M. A. Shay, and J. F. Drake, Phys. 
Rev. Lett. 95, 235002 (2005). 

16 D. A. Uzdensky, Astrophys. J. 671, 2139 (2007). 

17 P. A. Sweet, in Electromagnetic Phenomena in Cos- 
mical Physics, edited by B. Lehnert, IAU Symposium 
(Cambridge University Press, 1958), Vol. 6, p. 123. 

18 E. N. Parker, J. Geophys. Res. 62, 509 (1957). 

19 B. N. Rogers, R. E. Denton, J. F. Drake, and M. A. 
Shay, Phys. Rev. Lett. 87, 195004 (2001). 

20 N. F. Loureiro, A. A. Schekochihin, and S. C. Cowley, 
Phys. Plasmas 14, 100703 (2007). 

21 R. Samtaney, N. F. Loureiro, D. A. Uzdensky, A. A. 
Schekochihin, and S. C. Cowley, Phys. Rev. Lett. 
103, 105004 (2009). 



22 A. Bhattacharjee, Y. Huang, H. Yang, and 
B. Rogers, Phys. Plasmas 16, 112102 (2009). 

23 P. A. Cassak, M. A. Shay, and J. F. Drake, Phys. 
Plasmas 16, 120702 (2009). 

24 Y. Huang and A. Bhattacharjee, Phys. Plasmas 17, 
062104 (2010). 

25 D. A. Uzdensky, N. F. Loureiro, and A. A. Schekochi- 
hin, Phys. Rev. Lett. 105, 235002 (2010). 

26 W. H. Matthaeus and S. L. Lamkin, Phys. Fluids 29, 
2513 (1986). 

27 G. Kowal, A. Lazarian, E. T. Vishniac, and 
K. Otmianowska-Mazur, Astrophys. J. 700, 63 
(2009). 

28 G. Lapenta, Phys. Rev. Lett. 100, 235001 (2008). 

29 N. F. Loureiro, D. A. Uzdensky, A. A. Schekochihin, 
S. C. Cowley, and T. A. Yousef, Mon. Not. R. Astron. 
Soc. 399, L146 (2009). 

30 J. Birn, J. F. Drake, M. A. Shay, B. N. Rogers, 
R. E. Denton, M. Hesse, M. Kuznetsova, Z. W. Ma, 
A. Bhattacharjee, A. Otto, et al., J. Geophys. Res. 
106, 3715 (2001). 

31 P. Ricci, J. U. Brackbill, W. Daughton, and 

G. Lapenta, Phys. Plasmas 11, 4102 (2004). 

32 W. Daughton, J. Scudder, and H. Karimabadi, Phys. 
Plasmas 13, 072101 (2006). 

33 J. F. Drake, M. A. Shay, and M. Swisdak, Phys. 
Plasmas 15, 042306 (2008). 

34 W. Daughton, V. Roytershteyn, B. J. Albright, 

H. Karimabadi, L. Yin, and K. J. Bowers, Phys. Rev. 
Lett. 103, 065004 (2009). 

35 W. Daughton, V. Roytershteyn, B. J. Albright, 
H. Karimabadi, L. Yin, and K. J. Bowers, Phys. Plas- 
mas 16, 072117 (2009). 

36 R. D. Hazeltine, C. T. Hsu, and P. J. Morrison, Phys. 
Fluids 30, 3204 (1987). 

37 T. J. Schep, F. Pegoraro, and B. N. Kuvshinov, Phys. 
Plasmas 1, 2843 (1994). 

38 R. Fitzpatrick and F. Porcelli, Phys. Plasmas 11, 
4713 (2004). 

39 R. Fitzpatrick and F. Porcelli, Phys. Plasmas 14, 
049902 (2007). 

40 F. L. Waelbroeck, R. D. Hazeltine, and P. J. Morri- 
son, Phys. Plasmas 16, 032109 (2009). 

41 E. Tassi, P. J. Morrison, D. Grasso, and F. Pegoraro, 
Nucl. Fusion 50, 034007 (2010). 

42 W. Wan, Y. Chen, and S. E. Parker, Phys. Plasmas 
12, 012311 (2005). 



27 



43 B. N. Rogers, S. Kobayashi, P. Ricci, W. Dorland, 
J. Drake, and T. Tatsuno, Phys. Plasmas 14, 092110 
(2007). 

44 X. Y. Wang, Y. Lin, L. Chen, and Z. Lin, Phys. 
Plasmas 15, 072103 (2008). 

45 A. Perona, L.-G. Eriksson, and D. Grasso, Phys. 
Plasmas 17, 042104 (2010). 

46 R. Numata, W. Dorland, G. G. Howes, 
N. F. Loureiro, B. N. Rogers, and T. Tatsuno, 
e-print larXiv: 1107.58421 

47 E. A. Frieman and L. Chen, Phys. Fluids 25, 502 
(1982). 

48 A. A. Schekochihin, S. C. Cowley, W. Dorland, G. W. 
Hammett, G. G. Howes, G. G. Plunk, E. Quataert, 
and T. Tatsuno, Plasma Phys. Control. Fusion 50, 
124024 (2008). 

49 A. A. Schekochihin, S. C. Cowley, W. Dorland, G. W. 
Hammett, G. G. Howes, E. Quataert, and T. Tat- 
suno, Astrophys. J. Suppl. 182, 310 (2009). 

50 E. Cafaro, D. Grasso, F. Pegoraro, F. Porcelli, and 
A. Saluzzi, Phys. Rev. Lett. 80, 4430 (1998). 

D. Grasso, F. Califano, F. Pegoraro, and F. Porcelli, 
Phys. Rev. Lett. 86, 5051 (2001). 

52 K. Bowers and H. Li, Phys. Rev. Lett. 98, 035002 
(2007). 

53 H. J. de Blank, Phys. Plasmas 8, 3927 (2001). 

54 H. J. de Blank and G. Valori, Plasma Phys. Control. 
Fusion 45, A309 (2003). 

55 T. V. Liseikina, F. Pegoraro, and E. Y. Echkina, 
Phys. Plasmas 11, 3535 (2004). 

56 F. Pegoraro, T. Liseikina, and E. Yu Echkina, 
Transp. Th. Stat. Phys. 34, 243 (2005). 

57 M. Yamada, Y. Ren, H. Ji, J. Breslau, S. Gerhardt, 
R. Kulsrud, and A. Kuritsyn, Phys. Plasmas 13, 
052119 (2006). 

58 A. Zocco, L. Chacon, and A. N. Simakov, Phys. Plas- 
mas 16, 110703 (2009). 

59 T. Tatsuno, W. Dorland, A. A. Schekochihin, G. G. 
Plunk, M. Barnes, S. C. Cowley, and G. G. Howes, 
Phys. Rev. Lett. 103, 015003 (2009). 

60 S. I. Braginskii, in Reviews of Plasma Physics, edited 
by M. A. Leontovich (Consultants Bureau, New 
York, 1965), vol. 1, pp. 205-311. 

61 G. G. Howes, S. C. Cowley, W. Dorland, G. W. Ham- 
mett, E. Quataert, and A. A. Schekochihin, Astro- 
phys. J. 651, 590 (2006). 

62 H. P. Furth, J. Killeen, and M. N. Rosenbluth, Phys. 
Fluids 6, 459 (1963). 

63 B. Coppi, R. Galvao, R. Pellat, and M. N. Rosen- 
bluth, Sov. J. Plasma Phys. 11, 226 (1975). 

64 G. Ara, B. Basu, B. Coppi, G. Laval, M. N. Rosen- 
bluth, and B. V. Wadell, Ann. Phys. 112, 443 (1978). 

65 J. F. Drake, Phys. Fluids 21, 1777 (1978). 

66 J. F. Drake and Y. C. Lee, Phys. Fluids 20, 1341 
(1977). 

67 S. C. Cowley, R. M. Kulsrud, and T. S. Hahm, Phys. 
Fluids 29, 3230 (1986). 

68 F. Pegoraro and T. J. Schep, Plasma Phys. Control. 
Fusion 28, 647 (1986). 

69 F. Pegoraro, F. Porcelli, and T. J. Schep, Phys. Flu- 
ids B 1, 364 (1989). 



70 A. Y. Aydemir, Phys. Fluids B 3, 3025 (1991). 

71 G. Laval, R. Pellat, and M. Vuillemin, Plasma Phys. 
Control. Nucl. Fusion Res. 2, 259 (1965). 

72 B. Basu and B. Coppi, Phys. Fluids 24, 465 (1981). 

73 F. Porcelli, Phys. Rev. Lett. 66, 425 (1991). 

74 L. Zakharov and B. Rogers, Phys. Fluids B 4, 3285 
(1992). 

75 B. N. Kuvshinov, Plasma Phys. Control. Fusion 36, 
867 (1994). 

76 V. V. Mirnov, C. C. Hegna, and S. C. Prager, Phys. 
Plasmas 11, 4468 (2004). 

77 V. M. Vasyliunas, Rev. Geophys. 13, 303 (1975). 

78 J. A. Wesson, Nucl. Fusion 30, 2545 (1990). 

79 S. V. Bulanov, F. Pegoraro, and A. S. Sakharov, 
Phys. Fluids B 4, 2499 (1992). 

80 M. Ottaviani and F. Porcelli, Phys. Rev. Lett. 71, 
3802 (1993). 

81 W. Gekelman, H. Pfister, Z. Lucky, J. Bamber, 

D. Leneman, and J. Maggs, Rev. Sci. Instr. 62, 2875 
(1991). 

82 G. Saibene, N. Oyama, J. Lonnroth, Y. Andrew, 

E. de la Luna, C. Giroud, G. T. A. Huysmans, Y. Ka- 
mada, M. A. H. Kempenaars, A. Loarte, et al., Nucl. 
Fusion 47, 969 (2007). 

83 L. Spitzer and R. Harm, Phys. Rev. 89, 977 (1953). 

84 P. Helander and D. J. Sigmar, Collisional Trans- 
port in Magnetized Plasmas (Cambridge University 
Press, Cambridge, 2002). 

85 S. C. Cowley, PhD Thesis, Princeton University, 
1985. 

86 J. A. Krommes, Phys. Rep. 360, 1 (2002). 

87 A. Lenard and I. B. Bernstein, Phys. Rev. 112, 1456 
(1958). 

88 I. G. Abel, M. Barnes, S. C. Cowley, W. Dorland, 
and A. A. Schekochihin, Phys. Plasmas 15, 122509 
(2008). 

89 H. R. Strauss, Phys. Fluids 19, 134 (1976). 

90 I. G. Abel, G. G. Plunk, E. Wang, M. Barnes, 
S. C. Cowley, W. Dorland, and A. A. Schekochihin, 
Plasma Phys. Control. Fusion, submitted (2011). 

91 F. C. Grant and M. R. Feix, Phys. Fluids 10, 696 
(1967). 

92 T. P. Armstrong, Phys. Fluids 10, 1269 (1967). 

93 G. W. Hammett, M. A. Beer, W. Dorland, S. C. 
Cowley, and S. A. Smith, Plasma Phys. Control. Fu- 
sion 35, 973 (1993). 

94 S. E. Parker and D. Carati, Phys. Rev. Lett. 75, 441 
(1995). 

95 H. Sugama, T.-H. Watanabe, and W. Horton, Phys. 
Plasmas 8, 2617 (2001). 

96 G. G. Howes, W. Dorland, S. C. Cowley, G. W. Ham- 
mett, E. Quataert, A. A. Schekochihin, and T. Tat- 
suno, Phys. Rev. Lett. 100, 065004 (2008). 

97 N. F. Loureiro, A. Zocco, and A. A. Schekochihin, in 
preparation (2011). 

98 A. Kolmogorov, Dokl. Akad. Nauk SSSR 30, 301 
(1941). 

99 J. F. Drake and R. G. Kleva, Phys. Rev. Lett. 66, 
1458 (1991). 

100 G. Valori, D. Grasso, and H. J. de Blank, Phys. Plas- 
mas 7, 178 (2000). 



28 



D. Borgogno, D. Grasso, F. Porcelli, F. Califano, 

F. Pegoraro, and D. Farina, Phys. Plasmas 12, 
032309 (2005). 

G. G. Plunk, S. C. Cowley, A. A. Schekochihin, and 
T. Tatsuno, J. Fluid Mech. 664, 407 (2010). 

B. D. Fried and S. D. Conte, The Plasma Dispersion 

Function (Academic Press, New York, 1961). 

L. Chacon, A. N. Simakov, V. S. Lukin, and 

A. Zocco, Phys. Rev. Lett. 101, 025003 (2008). 

T. M. Antonsen and B. Coppi, Phys. Lett. A 81, 335 

(1981). 

J. F. Drake, J. T. M. Antonsen, A. B. Hassam, and 
N. T. Gladd, Phys. Fluids 26, 2509 (1983). 
T.-H. Watanabe and H. Sugama, Phys. Plasmas 11, 
1476 (2004). 

It is probably also essential, for a similar reason, that 
the model be three-dimensional; restriction to exact 
two-dimensionality imposes many additional conser- 
vation laws (see Appendix |A")) . 

It can also be thought to describe the nonlinear sit- 
uation, rather than the linearized one, if Eq. (|95H is 



first transformed into the Lagragian frame moving 
with the plasma and the parallel coordinate mea- 
sured along the exact field line. This, however, is not 
an essential point as the phase mixing mechanism we 
are considering is fundamentally linear. 

The spectrum with the sa me m dependence as (I105|) 
was first derived in Ref. 1 1071 for electrostatic ITG 
turbulence in tokamaks, although, unlike those au- 
thors, we have not invoked any nonlinear dynamics 
in its derivation. 

1 Note, however, the claim in Ref. [S3 that the inclusion 
of electron kinetics alleviates somewhat the problem 
of formation of singularities in their 2D Hamiltonian 
fluid-kinetic model. 

2 This reduction is accomplished by setting the elec- 
tron density gradient (1/a) and, therefore, cj, = 
kycTo e /eBoa to zero in all of their equations ex- 
cept Eq. (25), where (w*/cu) 2 //3 P = d 2 /2S 2 , because 
the density gradient in a;* cancels with the one in 
P = (4im 0e T 0e /B 2 )L 2 a /a 2 . 



